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1.0  INTRODUCTION 


High  Power  Microwave  (HPM)  research  of  the  past  10  to  15  years  has  focused  on 
developing  high  peak  power  (  >  1  GW),  high  energy  (  >  1  kJ)  sources  which  have 
moderate  frequency  tunability  (Refs.  1  to  20).  Repetitive  pulse  power 
operation  Ref.  13)  is  also  a  desired  parameter,  however  this  report  will  not 
address  the  rep-rate  issue.  The  standard  research  approach  is  to  develop  a 
concept,  model  the  instability  on  a  computer,  and  then  build  and  rebuild  the 
device  until  the  desired  performance  is  obtained.  This  procedure  typically 
requires  rebuilding  the  experimental  device  several  times,  which  is  typically 
expensive.  The  research  has  focused  on  how  this  process  may  be  improved. 
Specifically,  by  working  with  the  various  calculational  tools  (MAGIC,  SOS,  , 
EGUN,  etc.)  what  important  aspects  need  to  be  incorporated  into  these  tools  to 
enable  the  researcher  to  produce  better  designs  prior  to  "cutting  metal". 

The  technique  has  been  to  take  an  existing,  demonstrated  HPM  device  and 
compare  calculation  and  experiment  to  find  what  physics  may  not  have  been 
included  in  the  simulation,  but  may  be  essential  for  agreement  between  theory 
and  experiment.  Obviously,  these  codes  have  been  benchmarked  and  their 
validity  established,  however  the  rigor  of  the  solution  has  been  limited  to 
establishing  trends  but  not  numerical  design  accuracy.  Numerical  design 
accuracy  is  clearly  required  to  make  these  codes  true  source  design  tools. 

Based  on  the  published  results  of  the  Relativistic  Klystron  Amplifier  (RKA) 
(Ref.  4)  and  the  future  ability  to  phase  lock  amplifiers  for  HPM  arrays  we 
chose  the  RKA  as  a  prototype  HPM  system  to  begin  a  study  of  HPM  amplifiers. 
This  report  will  document  the  experimental  progress  of  the  Annular  Beam 
Amplifier  (ABA)  during  the  past  18  months,  and  the  success  simulating  the 
amplifier  (Fig.  1).  The  discussion  will  include  the  concept  for  building  a 
HPM  amplifier,  design  of  the  foilless  diode,  demonstration  of  the  electron 
beam  modulation,  and  the  insertion  of  assorted  diagnostics  to  compare  with  the 
computer  codes.  The  report  will  conclude  with  a  discussion  of  the  upcoming 
experiments.  For  future  reference  sample  input  files  for  the  various  computer 
codes  are  included  as  appendices. 
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2.0  EXPERIKENTAL  APPARATUS 


The  choice  of  the  RKA  allowed  one  to  break  the  problem  into  pieces  which  could 
be  individually  studied,  e.g.  the  vacuum  diode  to  provide  the  electron  beam 
power  to  be  converted  into  the  microwave  signal.  We  follow  the  procedure  of 
Reference  2  to  validate  coupling  the  signal  into  the  modulation  cavity.  One 
then  observes  the  growth  of  the  radiofrequency  (rf)  current  as  the  modulated 
beam  propagates  from  the  modulation  cavity.  At  the  peak  of  the  rf  current  on 
the  electron  beam  one  places  the  first  idler  cavity.  Then  again  the  peak  of 
the  rf  current  is  determined  and  another  idler  cavity  to  provide  further 
modulation,  or  an  extractor  is  placed  to  separate  the  rf  wave  from  the  beam. 

A  number  of  diagnostics  are  required  for  comparison  with  the  codes,  e.g. 
cavity  electric  field  probes  to  monitor  the  gap  voltages,  current  diagnostics 
to  observe  the  beam  optics  and  beam  losses,  and  beam  collector  to  observe  the 
cross-sectional  area  of  the  electron  beam. 

2.1  DESIGN  OF  AMPLIFIER 

To  build  flexibility  into  the  amplifier,  the  microwave  circuit  was  designed  to 
be  modular  and  to  be  insertable  within  a  vacuum  vessel.  This  allows  for 
future  investigations  with  other  microwave  circuits  without  needing  to  rebuild 
the  vacuum  vessel.  Tolerances  were  quite  restrictive  (0.001")  in  machining 
the  various  stainless  steel  parts.  Copper  finger  stock  was  used  to  make  good 
rf  electrical  contact.  This  was  done  to  maximize  the  quality  factor  (Q) .  As 
stated  above,  the  majority  of  the  parts  were  machined  from  stainless  steel  to 
allow  diffusion  of  a  pulsed  magnetic  field  (20  ms  rise  time)  through  the 
vacuum  vessel /microwave  circuit.  The  design  of  the  magnetic  field  coils  will 
be  discussed  below. 

The  guiding  philosophy  of  this  effort  has  been  to  make  all  components  as 
flexible  as  possible.  This  was  done  for  two  reasons.  First,  a  typical  set  of 
parameters  for  research  devices  do  not  match  requirements  of  other  users. 
Secondly,  the  computer  simulations  might  show  a  radical  design  change  to  be 
necessary.  This  was  found  to  be  true  early  in  the  research  program  when  the 
origi.ial  L.h..c..i  quarter  vv;..vo length  (K?.)  mode  modulating  cavity  supported  self¬ 
oscillation  in  the  quarter  wavelength  (ViX)  mode.  The  solution  was  to  complete 
the  simulation  with  ViK  cavities,  with  the  plan  to  suppress  the  fundamental 
mode  in  the  experiment  if  required.  This  lower  frequency  operation  has  been 
observed  by  Friedman  (Ref.  2) . 

The  cathode  placement  should  allow  generation  of  a  large  electron  beam  current 
pulse,  constrained  to  be  less  than  the  Child-Langmuir ,  space-charge  limiting 
current  in  the  beam  line.  This  removes  the  potential  of  torming  a  virtual 
cathode  (Ref.  19).  The  requirement  for  high  beam  power  requires  the  beam  to 
travel  close  to  the  beam  drift  pipe  wall.  The  beam  trajectory  limits  the 
maximum  tolerable  variation  of  the  axial  and  radial  magnetic  field.  Also, 
minimal  electron  beam  loss  through  the  microwave  circuit  is  required  to 
maximize  the  rf  power  and  to  reduce  plasma  formation  within  the  microwave 
circuit.  All  the  above  requirements  translate  to  the  highest  power/energy 
electron  beam  for  conversion  to  a  HPM  pulse. 

The  initial  solution  was  a  aluminum  "cookie  cutter"  cathode  of  outer  radius 
7.0 -cm  in  a  7.65-cm  beam  line  immersed  in  a  uniform  DC  axial  magnetic  field. 
The  second  cathode  (Fig.  2)  was  a  stainless  steel  right  circular  truncated 
cone  immersed  in  a  pulsed  uniform  magnetic  field,  with  electron  emission  from 
the  edge  of  the  base  of  radius  7.0  cm.  Both  cathodes  suffered  from  equal 
emitted  currents  parallel  and  anti-parallel  to  the  symmetry  axis.  This 
resulted  in  only  half  of  the  electron  beam  power  going  into  the  microwave 
amplifier  circuit. 
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Figure  2  Drawing  of  the  amplifier  using  the  inverted  truncated  cone  cathode. 
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The  present  cathode  is  being  tested  (Fig.  3).  The  focusing  magnetic  field 
causes  the  electron  beam  to  be  focused  to  a  smaller  radius  than  the  emission 
radius.  The  larger  diameter  of  the  cathode  provides  a  collector  for  electrons 
that  are  emitted  away  from  the  amplifier  cavities.  Since  the  electrons  are 
collected  at  the  cathode  potential  there  are  no  resultant  particles  created 
which  may  cause  problems  of  secondary  electrons  and  plasma  formation,  as  one 
would  expect  when  the  electrons  are  collected  at  ground  potential.  The 
cathode  has  been  successfully  simulated  using  EGUN  and  MAGIC  when  the  proper 
B.  and  B,  fields  are  used. 

The  microwave  circuit  and  drift  spaces  were  designed  to  allow  positional 
flexibility  during  the  project.  All  internal  sections  are  machined  to  0.001 
inch  and  finger  stock  is  used  to  provide  good  electrical  contact  (quality 
factors  of  700  to  800  were  easily  obtained).  All  sections  slide  within  the 
vacuum  vessel.  This  allows  easy  adjustments  in  cavity  dimensions,  gap  spaces, 
and  drift  spaces.  These  parameters  must  be  optimized  to  allow  the  highest 
efficiency  of  beam  power  to  rf  power  conversion. 

Not  shown  in  Figure  3  is  the  double  stub  tuner  built  to  provide  impedance 
matching  between  the  external  oscillator  and  the  B-dot  loop  used  to  excite  the 
modulation  cavity.  The  selection  of  a  B-dot  loop  with  a  vacuum  co-axial 
transmission  line  was  necessitated  by  the  space  available  to  inject  power  into 
the  amplifier.  Previous  experiments  have  used  waveguide  feeds  with  various 
iris  couplers  (Refs.  4,10,13).  The  iris  provides  the  same  impedance  matching 
as  obtained  with  the  double  stub  tuner.  A  double  stub  tuner  was  required 
since  proper  placement  of  a  single  stub  tuner  relative  to  the  B-dot  coupling 
loop  was  not  possible.  Calibration  data  will  be  presented  in  Section  4. 

Fazio  (Ref.  13)  reported  that  overcoupling  of  the  rf  oscillator  to  the 
modulator  cavity  was  required  due  to  electron  beam  loading  of  the  cavity. 

This  phenomena  has  been  observed  in  simulation  and  experiment  and  will  be 
discussed  later.  Data  are  presented  in  Section  4  which  show  how  this 
overcoupling  is  achieved  with  the  double  stub  tuner. 

The  drift  space  containing  an  axial  and  azimuthal  array  of  B-dot  probes,  shown 
in  Figure  3  following  the  idler  cavity,  is  used  to  measure  the  modulation 
generated  by  both  cavities.  An  array  of  20  zero  area  B-dot  loops  and  a 
Rogowskii  loop  are  used  to  monitor  beam  propagation  and  beam  modulation.  The 
array  is  composed  of  4  axial  lines  of  5  B-dots  each.  As  shown  in  Fig. 3  the 
lines  are  separated  in  angle  by  90°,  and  two  lines  are  axially  matched  while 
the  other  two  lines  are  positioned  to  be  centered  between  the  first  two  lines. 
This  provides  determination  of  the  azimuthal  symmetry  and  1-cm  resolution  on 
the  beam  modulation.  With  the  idler  cavity  removed  we  simply  slide  the 
"diagnostic  package"  toward  the  cathode,  causing  the  left  wall  of  the 
"diagnostic  package"  to  form  the  final  wall  of  the  modulating  cavity.  By 
using  the  same  probes  greatly  eases  the  calibration  requirements  and  data 
reduct  ion . 

An  rf  extractor  design  for  this  experiment  has  not  been  finalized.  Fazio 
(Ref.  13)  has  been  working  on  an  extractor,  however  the  Q  ot  the  extractor  at 
a  value  <  10  allows  one  to  develop  sufficient  voltage  to  reflect  the  electron 
beam  back  to  the  idler  cavity,  or  to  overvolt  the  extractor  gap  causing  an  rf 
arc.  At  this  point  the  amplifier  suffers  from  breakdown  and  the  rf  pulse 
terminates.  The  primaiy  problem  with  the  extractor  is  the  plasma  and 
secondary  particles  created  by  collecting  the  energy  of  the  unmodulated 
portion  of  the  electron  beam.  This  problem  is  not  restricted  to  HPM 
amplifiers;  any  HPM  device  which  requires  a  long  electron  beam  pulse,  long 
being  defined  as  :•  500  ns. 
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2.2  DIAGNOSTICS  AND  CALIBRATIONS 


In  completing  this  experiment,  we  have  tried  to  use  sufficient  diagnostics  to 
validate  critical  data  and  to  allow  comparison  with  simulation.  This  has 
resulted  in  using  E-dot  probes  in  the  rf  cavities  to  monitor  the  electric 
field  in  the  cavities  and  modulating  gaps,  B-dot  probes  in  the  drift  spaces  to 
monitor  the  beam  current  modulation  as  a  measure  of  the  amplifier  beam  power 
to  rf  power  efficiency,  and  a  Rogowskii  current  monitor  for  the  current 
exiting  the  drift  space. 

The  E-dot  probes  (Ref.  21)  are  cut  flush  tc  the  wall,  and  made  from  SMA 
bulkhead  feedthroughs.  These  probes  are  calibrated  with  a  ruler,  and  also 
cros.'  calibrated  with  the  oscillator  input  and  reflected  power  monitors  to 
validate  that  all  the  oscillator  powi  r  is  injected  in  the  cavity  during  the 
electron  beam  pulse.  This  coupling  coefficient  is  the  Sji  parameter  found  on 
a  vector  network  analyzer.  Those  E-dot  probes  are  also  used  to  determine  the 
axial  electric  f ield'potential  applied  to  the  modulator  gap  during  the 
presence  of  the  electron  beam.  The  ratio  of  the  radial  electric  field  at  the 
probe  to  the  axial  gap  electric  field  at  the  radius  of  the  beam  may  be 
estimated  by  SHY  (Ref.  22)  and  MAGIC  (Ref.  23).  The  accuracy  of  this  estimate 
will  have  to  be  determined  in  future  measurements  of  the  beam  modulation  or 
extracted  power. 

The  B-dot  probes  are  zero  area  probes  based  on  the  design  of  Voss'.  The  B- 
dot  probes  are  calibrated  on  the  Transmission  Line  Calibration  Fixture  ("Tin 
Man*)  a  50  Q  air  coaxial  transmission  line.  Tin  Man  allows  determination  of 
the  effective  area  of  the  B-dot  probe  (see  Appendix  A),  which  is  ideally  found 
to  be  independent  of  the  frequency  of  the  oscillator  up  to  some  maximum 
determined  by  the  geometry  of  the  probe.  This  area  may  also  be  determined  by 
pulsing  currents  on  a  rod  on  the  geometric  axis  of  the  “diagnostic  package*  or 
by  generating  an  electron  beam  and  propagating  the  beam  past  the  B-dot  probes 
and  the  exit  Rogowskii.  Without  an  external  oscillator  signal  injected  into 
the  modulation  cavity  these  diagnostics  should  agree,  modulo  the  geometric 
centering  of  the  electron  beam. 

The  anticipated  results  are  that  the  3  different  calibration  techniques  should 
provide  the  same  effective  area  for  each  B-dot  probe.  The  results  from  “Tin 
Man*  were  a  factor  of  1/V2  lower  than  found  from  beam  shots  or  short  shots. 
This  was  due  to  an  implicit  assumption  that  “Tin  Man*  uses  CW,  or  RMS  power 
rather  than  the  pulse  value  used  in  beam  or  short  shots.  Graphs  of  the  data 
with  and  without  the  factor  of  1/V2  show  this  comparison  (Fig.  4). 

A  primary  consideration  of  the  experimental  design  was  that  enough  diagnostics 
be  available  to  monitor  the  important  physics  of  the  amplifier.  This  includes 
the  requirement  that  the  diagnostics  allow  for  comparison  with  simulation. 
Also, the  diagnostic  techniques  must  be  self-consistent.  That  is,  if  we 
observe  a  certain  potential  in  the  gap  with  the  modulation  cavity  E-dot  probe, 
then  the  induced  beam  modulation  found  with  the  “diagnostic  package*  B-dot 
array  should  be  consistent  with  simulation  results. 


*  D.E.  Voss,  private  communication 
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Ab-rms  ^  As-rms  ♦  A-tInman 


Figure  4  Comparison  of  B-dot  probe  effective  areas. 
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3.0  SIMUIiATIONS 


An  assortment  of  simulations  or  calculations  have  been  completed  in  preparing 
for  the  design  of  this  experiment.  Sample  results  of  these  simulations  are 
presented  in  this  section.  Complete  input  data  files  and  additional  input  are 
listed  in  the  appendices. 

3.1  SUPERFZSH 

The  POISSON/  (Ref.  22)  group  of  codes  were  used  to  simulate  the  “cold"  steady- 
state  response  of  the  amplifier  modulation  and  idler  cavities.  The 
calculation  provides  information  on  the  resonant  modes  of  the  microwave 
cavities  and  measurements  of  electric  field  amplitudes  throughout  the  cavity. 
The  code  group  is  very  well  documented  in  the  users  manuals,  and  the  reader  is 
referred  to  them. 

To  use  the  group  of  codes,  one  begins  with  a  geometry  of  the  device.  This 
input  is  first  processed  by  the  code  AUTOMESH.  Two  output  files  are  created. 
One  file  is  the  output  on  how  a  triangular  mesh  is  fit  to  the  problem  and  the 
second  is  an  input  file  for  the  code  LATTICE.  During  execution  of  LATTICE, 
the  AUTOMESH  output  file  is  read  and  the  user  is  prompted  for  input  for  the 
“CON"  array,  which  may  be  used  to  change  various  defaults  such  as  the  problem 
boundary  conditions,  problem  symmetry,  etc,  LATTICE  also  produces  two  output 
files,  the  first  being  an  ASCII  file  summarizing  the  execution  of  the  program 
LATTICE,  and  the  second  is  an  input  file  for  the  codes  PSFPLOT  and  .  PSFPLOT 
may  be  used  to  verify  the  geometry,  mesh  relaxation,  and  to  plot  the  electric 
field  profile  of  the  cavity.  is  executed  to  determine  the  resonant  frequency 
of  the  cavity  mode.  An  initial  guess  for  the  resonant  frequency  is  used  to 
select  which  mode  the  code  iterates  toward.  The  accuracy  of  the  guess 
determines  how  hard  the  code  must  work  to  find  a  solution.  One  may  also  have 
search  a  range  of  frequencies  to  look  for  the  number  of  modes  supported  by  the 
cavity  within  that  frequency  range.  A  postprocessor  code  SHY  may  be  used  to 
determine  ratios  of  electric  field  components  throughout  the  geometry.  This 
was  used  in  conjunction  with  electric  field  probes  to  determine  the  modulation 
voltage  established  by  the  modulation  cavity.  Sample  input  is  listed  in 
Appendix  B.  Results  for  the  lowest  two  modes  of  the  ABA  are  shown  in  Figure 
5. 

3.2  MAGNET 

The  magnetic  field  coils  were  initially  designed  using  the  MAGNET*  code 
developed  by  John  Freeman  of  Sandia  National  Laboratory.  This  code  requires  a 
definition  of  the  geometry  of  the  magnetic  field  coils,  specifically  the  left 
and  right  boundaries  of  each  coil,  the  number  of  turns  per  layer  and  number  of 
layers  for  each  coil,  and  the  current  flowing  in  each  layer  (see  sample  input 

deck  in  Appendix  C.).  By  setting  the  current  to  1  Ampere  the  code  results 

were  in  units  of  Gauss  per  Ampere.  The  code  calculates  the  spatial  variation 
of  the  steady  state  axial  and  radial  magnetic  fields  along  various  radial 

contours  as  a  function  of  the  axial  coordinate.  This  requires  one  to  verify 

that  magnetic  diffusion  is  not  an  issue  when  using  pulsed  coils. 

The  large  number  of  turns  (550)  used  in  the  magnetic  field  coil  design 
required  a  better  way  to  define  the  coil  geometry  to  MAGNET.  A  short  Fortran 
program  was  written  to  generate  the  input  deck  for  MAGNET.  This  allowed 
faster  changes  to  the  input  deck,  and  sped  the  design  process.  For  this 
application  minimal  axial  non-uniformity  at  the  beam  radius  was  required  to 
achieve  beam  transport.  A  check  on  the  axial  uniformity  was  accomplished  by 
calculating  the  radial  magnetic  field  component  at  the  beam  radius  and 


*  M.D.  Haworth,  private  communication 
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plotting  this  component  against  the  axial  position.  Additionally,  the  quality 
of  the  magnetic  field  system  is  indicated  by  comparing  the  axial  variation  of 
Bj  on  axis  and  at  the  beam  radius (Figure  6).  The  radial  and  axial  fields 
computed  by  MAGNET  were  reformatted  into  arrays  and  used  as  input  for  the 
applied  magnetic  field  used  in  the  MAGIC  PIC  simulations  of  the  RKA  (Figure 
7)  . 

A  comparison  of  the  magnetic  flux  of  the  designed  magnet  was  done  with  the 
code  POISSON  (Ref.  22).  POISSON  calculates  the  magnetic  flux  distribution  as 
functions  of  r  and  z.  POISSON  belongs  to  the  same  group  of  codes  as  .  One 
defines  the  geometry  similarly  to  SUPERFISH,  however,  where  SUPERFISH  uses 
cylindrical  z-r  on  an  x-y  grid  to  determine  cylindrical  electromagnetic  modes, 
POISSON  uses  cylindrical  r-z  on  an  x-y  grid  for  cylindrical  magnetic  flux 
plots.  The  magnetic  flux  plots  were  also  used  to  aid  in  the  placement  of  the 
new  cathode  (Fig.  8). 

An  additional  consideration  in  the  design  of  the  magnetic  field  coils  is 
whether  or  not  a  pulsed  current  source  is  to  be  used  to  generate  the  magnetic 
field.  The  eddy  currents  generated  by  a  pulsed  magnetic  field  passing  through 
the  material  used  to  assemble  the  vacuum  vessel  cause  a  time  lag  of  the 
magnetic  field  amplitude  within  the  electron  beam  line.  The  code  MAGDIF*  was 
used  to  simulate  the  magnetic  flux  penetration  through  the  stainless  steel 
vacuum  vessel.  The  results  showed  that  the  conductivity  of  stainless  steel 
was  not  sufficient  to  reduce  the  magnetic  flux  by  the  presence  of  eddy 
currents.  A  sequence  of  time  snapshots  of  the  magnetic  flux  generated  by 
MAGDIF  are  shown  in  Figure  9  and  the  steady  state  results  in  Figure  10  which 
agree  well  with  Figure  8. 

3.3  EGUN 

The  trajectory  code  EGUN  (Ref.  24)  is  used  to  design  the  vacuum  diode.  EGUN 
determines  the  steady-state  current  emitted  by  the  cathode,  and  the  final 
trajectories  of  the  emitted  electrons  subject  to  the  applied  potentials  and 
magnetic  fields.  A  sample  of  the  result  of  the  diode  used  in  a  uniform 
magnetic  field  is  shown  in  Figure  11.  The  input  file  is  listed  in  the 
Appendix  D.  The  converging  magnetic  field  was  included  later  to  verify  the 
diode  design. 

EGUN  is  described  fairly  completely  in  the  user's  manual.  However,  this 
computational  tool  is  only  easily  used  by  an  experienced  user,  and  the  input 
deck  is  somewhat  counter-intuitive  in  the  method  of  defining  conductors.  The 
version  we  had  access  to  did  not  include  all  the  present  interfaces  to 
commercial  CAD  software  and  this  probably  eases  use  by  novice  users.  The 
definition  of  the  magnetic  field  is  also  nonintuitive  from  the  point  of 
defining  magnetic  fields  for  annular  electron  beams.  The  magnetic  field  is 
defined  at  each  mesh  point  on  axis  and  then  computed  for  the  beam  radius. 

This  feature  led  to  confusion  in  understanding  the  magnetic  field  components 
computed  for  the  beam  radius.  It  required  some  time  to  understand  why  the 
magnetic  field  did  not  resemble  any  magnetic  field  profiles  calculated  by 
other  codes.  The  beam  emission  results  were  only  understood  following 
comparison  with  MAGIC. 


*  G.  Marklin,  private  communication 
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(Gauss) 


Figure  6  Comparison  of  axial  magnetic  field  on  aocis  and  at  the  beam  radius. 
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Axial  Magnetic  Field 
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Figure  8  Sample  magnetic  flux  contour  plot  from  POISSON. 


Figure  9  Magnetic  flux  snapshots  from  MAGDIF. 


Figure  10  Steady  state  result  from  MAGDIF. 


3.4  RKA 


The  one-dimensional  (1-D)  buncher  code  RKA  was  written  to  allow  rapid 
simulation  of  the  amplifier  gap  space  and  drift  spaces  is  based  on  the  theory 
of  Solymar  (Ref.  25).  This  code  propagates  the  electron  becim  through  a  single 
modulating  gap,  with  a  specified  voltage,  frequency,  and  gap  extent.  The  beam 
space  charge  is  an  input  parameter  which  is  calculated  following  the  formalism 
of  Reference  26.  The  user  must  provide  as  parameters  the  voltage  applied  at 
the  gap,  and  the  space  charge  reduction  factor  or  plasma  reduction  factor 
(PRF) .  One  postulates  an  applied  voltage  on  the  first  gap,  and  the  code 
provides  the  beam  modulation  through  the  drift 'Space  following  the  gap.  The 
output  data  file  may  then  be  used  as  input  to  any  subsequent  modulating  gaps 
(which  may  be  either  externally  or  self-excited) .  This  program  allows  for  any 
number  of  modulating  gaps  to  be  used  to  modulate  the  electron  beam  by  changing 
nrun  to  1,  and  renaming  the  output  data  file  to  be  the  input  data 
file  for  the  subsequent  simulation  (see  Appendix  E) . 

A  simulation  for  the  predicted  beam  modulation  and  optimum  phase  bunching  of 
the  electron  beam  kinetic  energy  for  the  final  experimental  configuration  and 
parameters  are  shown  below  (Figs.  12  and  13).  The  table  of  input  values 
yielded  the  graphical  results  of  two  RKA  simulations.  The  first  begins  at  the 
input  to  the  first  gap  and  concludes  at  the  z  coordinate  of  the  second  gap. 

The  output  of  the  first  simulation  was  used  as  input  for  the  second 
simulation.  The  second  gap  voltage  is  initially  chosen  to  give  the  desired 
maximum  modulation  coefficient.  If  this  value  is  unrealistic  in  amplitude, 
then  one  must  redesign  the  experiment  to  remove  this  problem. 


Figure  12  RKA  result  of  modulation  coefficient  versus  axial  drift  distance. 
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Figure  13  Sample  result  from  RKA  showing  kinetic  energy  phase  bunching  versus  drift  distance. 


That  this  code  is  very  versatile  and  allows  one  to  rapidly  try  any  sequence  of 
modulation  cavities,  drift  distances,  and  appl ied/ inauced  modulating  voltages. 
The  main  limitation  is  that  the  applied  voltage  a  is  that  voltage  actually  at 
the  position  of  the  becim,  not  the  voltage  developed  in  the  gap.  This  may  be  a 
difficult  value  to  estimate.  The  value  of  the  PRF  is  not  determined  self- 
consistently  as  the  beam  is  modulated.  These  issues  are  minor,  because 
excellent  agreement  with  MAGIC  has  been  obtained. 
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3.5  MAGIC 


The  particle-in-cell  (PIC)  code  MAGIC  has  been  used  to  model  the  electron  beam 
dynamics  of  an  intense  annular  electron  beam  propagated  through  an  RKA  system. 
MAGIC  is  a  2^-D,  fully  electromagnetic  and  relativistic,  self-consistent  PIC 
code  (Ref.  23) .  The  major  limitation  is  that  one  spatial  dimension  must  be 
assumed  ignorable,  or  that  any  motion  is  isotropic  in  that  dimension.  For 
this  project  this  limitation  is  avoided  by  the  azimuthal  symmetry  of  the 
device.  To  investigate  if  any  azimuthal  dynamics  occur,  one  would  have  to 
employ  a  full  3-D  code  such  as  SOS  or  ISIS. 

MAGIC  has  been  used  to  simulate  the  propagation  of  the  electron  beam  through 
the  RKA  without  an  applied  voltage  on  the  first  gap.  This  demonstrated  that 
the  ABA  simulation  did  not  self -oscillate .  The  simulation  was  an  injection 
run  with  a  uniform  axial  magnetic  field  of  8  kG.  The  injection  simulation 
allows  one  to  complete  the  simulation  much  quicker,  however  the  accuracy  may 
be  suspect  as  self-consistency  is  sacrificed.  The  simulations  typically  only 
match  the  first  60  ns  of  the  300  ns  pulse.  This  is  due  to  the  inherent 
problem  of  numerical  noise  generated  by  the  finite  differencing  employed  to 
solve  the  system  of  differential  equations.  This  problem  is  not  only  with 
MAGIC  but  with  all  finite  difference  technique  codes.  The  limitation  on 
problem  time  also  impacts  the  ability  to  simulate  other  physics,  such  as 
plasma  production  on  grounded  surfaces  or  self-ionization  of  background  gas. 
These  phenomena  typically  occur  on  a  much  longer  time  scale  than  one  can 
simulate  with  present  computer  systems. 


A  problem  occurred  in  trying  to  complete  an  emission  simulation  employing  the 
original  cathode.  The  technique  of  gridding  the  problem  space  causes  "stair¬ 
stepping*  on  slanted  lines.  This  causes  enhanced  emission  from  the  corners  of 
the  jagged  lines.  If  one  tries  to  increase  the  number  of  mesh  points  to  allow 
better  replication  of  the  “real"  system,  problems  arise  in  execution  speed 
because  of  the  number  of  mash  points.  These  issues  are  being  raised,  because 
they  limit  ones  ability  to  make  the  code  predictive  instead  of  a  qualitative 
guide . 

A  uniform  axial  magnetic  field  was  used  in  the  initial  simulations.  Clearly 
that  does  not  accurately  reflect  the  magnetic  field  generated  by  the  pulse 
field  coils.  MAGNET  was  used  to  generate  a  series  of  radial  contours  of  the 
axial  and  radial  magnetic  fields  as  functions  of  axial  coordinate.  These 
contours  were  then  put  in  matrix  form  and  read  into  MAGIC.  This  allows 
simulation  of  the  magnetic  focusing  used  in  the  new  "Friedman/Fazio*  cathode 
(Refs.  6,13).  The  curving  magnetic  flux  lines  are  crucial  to  demonstrating 
that  the  full  16  )cA  of  electron  beam  current  goes  into  the  amplifier  and  not 
away  from  the  region  of  interest.  The  matrices  for  and  B^  are  shown  below: 
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Bj  (Tesla)  = 


0.3863 

0.3854 

0.3825 

0.3772 

0.3689 

0.3566 

0.3393 

0.3158 

0.2855 

0.4353 

0.4347 

0.4328 

0.429 

0.4225 

0.4121 

0.3957 

0.3712 

0.3368 

0.4867 

0.4866 

0.4862 

0.485 

0.4821 

0.4759 

0.4635 

0.441 

0.4039 

0.5383 

0.5389 

0.5405 

0.5431 

0.5458 

0.5473 

0.5441 

0.5299 

0.4949 

0.5875 

0.5888 

0.5929 

0.5999 

0.6099 

0.6226 

0.6359 

0.6424 

0.6234 

0.6314 

0.6334 

0.6396 

0.6508 

0.6681 

0.6938 

0.7303 

0.7776 

0.8135 

0.6676 

0.6699 

0.6773 

0.6909 

0.713 

0.7479 

0.8046 

0.9031 

1.091 

0.6944 

0.6968 

0.7041 

0.7174 

0.739 

0.7734 

0.8298 

0.9306 

1 . 14 

0.7123 

0.7142 

0.72 

0.7304 

0.7464 

0.7702 

0.8047 

0.8513 

0.8903 

0.7228 

0.724 

0.7276 

0.7335 

0.7415 

0.7508 

0.7588 

0.7575 

0.7269 

0.7284 

0.7289 

0.7304 

0.7321 

0.7328 

0.73 

0.7188 

0.6909 

0.6344 

0.7318 

0.7319 

0.7318 

0.7309 

0.7273 

0.718 

0.6976 

0.6583 

0.5917 

0.7347 

0.7346 

0.7342 

0.7326 

0.7285 

0.7191 

0.6994 

0.6609 

0.5924 

0.7378 

0.7379 

0.7378 

0.7374 

0.7361 

0.7325 

0.7236 

0.7012 

0.6445 

0.7408 

0.7411 

0.7418 

0.7433 

0.7461 

0.7514 

0.7613 

0.7766 

0.777 

0.743 

0.7433 

0.7446 

0.7473 

0.7529 

0.7645 

0.7898 

0.8489 

1.011 

0.7438 

0.7441 

0.7452 

0.7476 

0.7525 

0.7626 

0.7844 

0.8331 

0.9409 

0.7434 

0.7435 

0.7438 

0.7445 

0.7457 

0.7477 

0.7507 

0.7511 

0.7241 

0.7426 

0.7424 

0.7419 

0 .7406 

0.7377 

0.7311 

0.7167 

0.6846 

0.6166 

0.7423 

0.742 

0.7411 

0.7388 

0.7339 

0.7234 

0.7017 

0.6589 

0.5822 

0.7429 

0.7428 

0.7421 

0.7405 

0.7368 

0.7289 

0.7118 

0.6757 

0.6042 

0.7442 

0.7443 

0.7444 

0.7446 

0.7448 

0.7448 

0.7433 

0.7343 

0.6931 

0.7453 

0.7456 

0.7466 

0.7487 

0.753 

0.7617 

0.78 

0.8186 

0.8854 

0.7454 

0.7458 

0.7471 

0.7499 

0.7557 

0.7677 

0.7941 

0.8571 

1.048 

0.744 

0.7443 

0.7452 

0.7472 

0.7509 

0.7581 

0.77.7  3 

0.7982 

0.824 

0.7415 

0.7416 

0.7418 

0.742 

0.7417 

0.7402 

0.735 

U.7188 

0 . 6686 

0.7386 

0.7386 

0.7383 

0.7372 

0.7338 

0.7256 

0.7077 

0.6713 

0.6034 

0.7357 

0.7358 

0.7353 

0.7349 

0.7315 

0.7222 

0.7018 

0.6622 

0.5945 

0.7323 

0.7328 

0.7341 

0.7356 

0.7358 

0.7321 

0.7193 

0.6892 

0.6304 

0.7268 

0.728 

0.7313 

0.7368 

0.7439 

0.7515 

0.7566 

0.7509 

0.7152 

0.7168 

0.7186 

0.7242 

0  ."’34 

0.7492 

0.7713 

0.8023 

0.8418 

0.8678 

0.6996 

0.7019 

0.709 

0.7221 

0.7432 

0.7767 

0.8312 

0.9269 

1  .111 

0.6735 

0.6758 

0.6833 

0.697 

0.7192 

0.7545 

0.8124 

0.915 

1.126 

0.6381 

0.6401 

0.6465 

0.658 

0.6761 

0.7031 

0.7424 

0.7957 

0.8433 

0.5947 

0.5961 

0.6005 

0.6079 

0.6188 

0.6329 

0.6486 

0.6587 

0.6432 

0.5457 

0.5464 

0.5483 

0.5513 

0.5548 

0.5573 

0.5556 

0.543 

0.5088 

0.494 

0.494 

0.4938 

0.493 

0.4906 

0.4851 

0.4734 

0.4514 

0.414 

0.4424 

0.4418 

0.44 

0.4364 

0.4303 

0.4201 

0.4039 

0.3794 

0.3444 

0.3929 

0.392 

0.3891 

0.3839 

0.3757 

0.3635 

0.3461 

0.3224 

0.2915 

0.3469 

0.3458 

0.3424 

0.3365 

0.3276 

0.315 

0.2982 

0.2766 

0.25 

0.3052 

0.304 

0.3005 

0.2945 

0.2857 

0.2737 

0.2582 

0.2392 

0.2167 

0.268 

0.2669 

0.2635 

0.2578 

0.2495 

0.2385 
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One  then  observes  the  beam  modulation  build  up  as  the  beam  propagates  through 
the  ABA,  and  the  induced  voltage  developed  on  the  idler  cavity.  MAGIC  employs 
a  command  called  "antenna"  which  specifies  a  current  density  on  a  wire  dipole 
model.  To  compare  experiment  and  simulation  would  require  calibrating  the 
amplitude  to  know  the  power  required  from  the  external  oscillator.  Instead, 
we  monitor  the  rf  mode,  the  voltage  on  the  modulation  gap  and  infer  the  seutie 
voltage  from  an  electric  field  probe  in  the  cavity.  This  indicates  whether 
the  input  cavity  is  driven  with  sufficient  power  to  modulate  the  electron 
beam.  Sample  plots  from  MAGIC  of  the  voltage  on  the  first  and  second  gap  are 
shown  (Figs.  14  and  15). 


The  drop  in  voltage  (Fig.  14)  is  due  to  beam  loading  of  the  modulating  gap. 
This  effect  has  been  seen  by  Fazio,  and  will  be  discussed  later.  The 
amplitude  of  the  voltage  on  the  second  gap  corresponds  well  with  the  assumed 
voltage  used  in  the  second  gap  of  the  RKA  simulations.  A  representative  plot 
(Fig.  16)  of  the  beam  modulation  based  on  simulations  with  MAGIC  and  RKA  is 
shown.  Note  the  qualitative  and  quantitative  agreement. 
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Figure  14  Time  variation  of  modulation  gap  voltage  from  MAGIC. 
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Figure  16  Comparison  of  MAGIC  and  RKA  results  of  axial  variation  of  beam  modulation. 


4.0  EXFERXHENTAI  RESULTS 


4.1  OSCILLATOR  COUPLING/BEAM  LOADING 

To  develop  the  modulating  voltage  across  the  modulating  gap,  a  signal  from  an 
external  oscillatoris  requied.  A  500-kW,  2-p.sec  magnetron,  which  is  tunable 
from  1.25  GHz  to  1.35  GHz  was  employed.  We  use  a  12-dB  isolator  on  the  output 
of  the  magnetron  to  reduce  the  reflected  signal  coupling  into  the  magnetron. 

A  dual  waveguide  loop  coupler  is  used  to  monitor  the  input  power  to  the  cavity 
and  the  reflected  power  from  the  B-dot  coupling  loop  in  the  modulating  cavity. 
These  monitors  are  used  both  in  cold  tests  and  experiments  with  the  electron 
beam.  The  difference  in  the  two  waveguide  power  monitors  should  time 
correlate  with  power  detected  in  the  cavity  by  the  E-dot  probe. 

In  principle  one  could  construct  a  B-dot  loop  in  the  cavity  with  the  proper 
cross  sectional  area  and  conductor  size  to  match  the  waveguide  impedance.  In 
practice,  this  is  very  difficult  because  the  mode  within  the  cavity  must  be 
known  to  determine  the  area  required  to  match  the  real  part  of  the  impedance, 
and  adjusting  the  size  of  the  conductor  used  in  forming  the  loop  to  cancel  the 
reactive  impedance  is  impractical.  Also,  this  impedance  match  must  be  done 
for  each  frequency  used  in  the  amplifier.  A  solution  was  to  use  a  double  stub 
tuner  which  could  easily  be  adjusted  outside  the  experimental  system 
regardless  of  the  frequency  of  operation.  We  chose  a  double  stub  tuner  over  a 
single  stub  tuner  simply  due  to  our  inability  to  place  the  tuner  at  a  maximum 
in  the  line  standing  wave  ratio  (SWR)  pattern. 

Figures  17  and  18  show  data  for  the  reflection  and  transmission 
characteristics  of  the  B-dot  coupling  loop  into  the  modulation  cavity.  The 
reflection  is  shown  by  the  parameter,  and  the  transmission  is  shown  by  the 
S21  parameter  from  the  network  analyzer.  The  calibration  procedure  used 
standard  50  Q.  RG-214  cable  to  connect  the  analyzer  to  the  experiment.  The 
device  under  test  included  a  type-N  to  WR-650  adapter,  a  section  of  WR-650 
waveguide  with  a  dielectric  plate  used  to  make  an  air-vacuum  interface, 
another  WR-650  to  vacuum  coaxial  line,  the  double  stub  tuner,  the  coupling 
loop  within  the  cavity,  and  the  E-dot  cavity  probe.  This  allows  the  observed 
reflection  to  be  the  composite  of  all  the  connections  in  the  transmission 
line.  The  transmission  coefficient  (Fig.  18)  included  all  the  effective  areas 
of  the  coupling  loop  and  the  E-dot  probe.  The  reflection  was  typically  viewed 
on  a  Smith  Chart,  where  a  match  to  input  was  shown  as  an  intersection  with  the 
origin  of  the  chart. 

While  it  is  trivial  to  use  the  double  stub  tuner  to  match  the  B-dot  coupling 
loop  to  the  50  fl  waveguide  impedance  at  the  cavity  frequency,  we  have 
observed,  as  has  Fazio  (Ref.  13),  that  matching  the  cold  test  cavity  frequency 
does  not  provide  optimum  matching  during  the  presence  of  the  electron  beam. 

In  fact,  one  must  view  the  cold  cavity  impedance  as  shunted  by  the  electron 
becim  impedance.  This  combined  load,  which  is  dominated  by  the  lower  electron 
beam  impedance,  must  be  matched  to  the  waveguide /B-dot  launcher.  This 
parallel  impedance  causes  a  slightly  different  matching  condition  for  the  B- 
dot  coupler  into  the  modulation  cavity.  Two  methods  may  be  employed  to 
compensate  for  this  beam  loading  of  the  cavity.  First,  the  double  stub  tuner 
may  be  properly  adjusted  for  this  new  optimum  frequency.  Secondly,  if  the 
cavity  has  a  lower  beam  loaded  Q,  adjust  the  frequency  of  the  magnetron  by  a 
few  MHz  until  the  double  stub  tuner  provides  a  match  to  the  loaded  cavity. 

The  reader  should  recall  that  the  modulating  cavity  voltage  observed  in  the 
MAGIC  simulations  also  shows  this  beam  loading.  The  voltage  drop,  which 
occurs  after  the  10-ns  delay  to  fill  the  cavity,  is  a  manifestation  of 
electron  beam  loading. 
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Figure  17  Results  of  modulation  cavity  cold  tests  showing  impedance  matching  of  waveguide  to 
cavity. 


satisfied. 


4.2  BEAM  PROPAGATION 

One  of  the  first  experiments  to  be  completed  in  this  type  of  project  is  the 
generation  and  propagation  of  an  intense  electron  beam  through  the  microwave 
circuit  and  extractor  region  of  the  amplifier.  To  generate  an  intense,  or 
high  power,  electron  beam  requires  that  the  product  of  the  beam  voltage  and 
current  be  large,  typically  on  the  order  of  a  few  GW  of  power.  This 
requirement  specifies  the  voltage  or  impedance  of  the  electron  beam.  A 
constraint  is  that  the  electron  beam  current  must  be  below  the  space  charge 
limiting  current  for  the  beam  cross  section  in  the  beamline  of  the  vacuum 
vessel.  For  this  experiment  with  a  7.65-cm  radius  beam  line,  and  the  electron 
beam  between  6.6  to  7.1  cm,  one  finds  that  the  space  charge  limiting  current 
is  22  kA  for  a  400-kV  electron  beam.  As  the  beam  modulation  grows  or  as  the 
beam  passes  by  the  modulation  and  idler  gaps,  the  space  charge  limiting 
current  condition  is  violated.  However,  the  transit  time  of  the  gaps  is  small 
compared  to  the  rf  wavelength,  therefore  the  potential  for  self-oscillation  or 
reflection  of  electrons  is  minimized. 

A  series  of  experiments  were  conducted  with  both  cathode  geometries  to 
demonstrate  the  propagation  of  the  electron  beeim  through  the  microwave 
circuit.  The  original  cathode  in  a  uniform  magnetic  field,  the  modulating 
cavity  and  the  diagnostic  package  with  an  inserted  graphite  witness  plate  were 
used  to  image  the  electron  beam  at  various  axial  positions.  The  A-K  gap  in 
this  configuration  is  simply  the  radial  separation  of  the  cathode  and  wall. 

The  beam  propagated  for  a  maximum  of  55  cm.  The  witness  plate  showed  a  time 
integrated  image  of  the  electron  beam  which  was  of  5-mm  uniform  thickness,  and 
the  image  inner  radius  matched  the  outer  radius  of  the  original  cathode 
(Fig. 19).  The  only  occasion  which  the  witness  plate  showed  a  thicker  image 
was  when  the  plate  was  a  few  centimeters  from  the  cathode.  The  thickening  of 
the  beam  pattern  was  due  to  the  psuedo-ground  formed  by  the  graphite  plate. 


Figure  19  Sample  of  graphite  witness  plate  data  for  uniform  field  propagation. 
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A  second  series  of  experiments  was  done  with  the  second  cathode  located  in 
both  a  converging  and  uniform  magnetic  field.  In  this  geometry  the  A-K  gap  is 
primarily  axial,  but  there  is  a  small  radial  separation.  The  beam  is 
propagated  through  the  modulating  gap,  idler  gap,  and  the  diagnostic  package. 
The  total  axial  distance  is  65  to  70  cm.  Recall  the  magnetic  field  coils  are 
only  deemed  uniform  for  about  50  cm.  These  experiments  dealt  with  varying  the 
position  of  the  cathode  in  the  portion  of  the  magnetic  field  where  the  field 
changes  from  mainly  B,  to  and  the  size  of  the  A-K  gap.  The  goal  was  to 
determine  the  proper  positions  to  get  the  16  kA  beam  through  the  beeunline,  and 
to  have  the  witness  plate  indicate  that  the  beam  has  the  desired  thickness  and 
radius  (Fig.  20) .  This  set  of  experiments  also  provided  a  check  on  the 
accuracy  of  the  B-dot  probe  array  calibrations,  and  the  azimuthal  uniformity 
of  the  electron  beam. 

The  position  of  the  cathode,  in  order  to  provide  the  proper  witness  plate 
pattern  was  about  6  cm  further  from  the  anode  than  anticipated.  This  also 
required  fabrication  of  an  anode  extension  of  4  cm  such  that  the  A-K  gap 
remained  at  about  2  cm.  Sample  data  are  shown  in  Figure  21  for  the  case  of  a 
3-cm  anode  extension,  3-cm  A-K  gap,  and  the  magnets  in  the  original  position. 
Notice  that  the  current  diagnostics  indicate  about  2  kA  of  current  loss, 
however  the  witness  plate  indicates  that  the  beam  is  several  rom  from  the  wall. 
Inspection  of  the  “cathode  end"  of  the  anode  found  several  burn  marks 
indicating  beam  scrape-off  at  the  anode.  The  B-dots  of  the  diagnostic  package 
agree  with  the  exit  Rogowskii  indicating  no  additional  loss  as  the  beam 
travels  through  the  structure. 
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positions. 


4.3  BEAM  MODUIiATION 

While  using  the  original  cathode  system,  there  was  a  brief  opportunity  to 
attempt  to  modulate  the  electron  beam.  The  data  are  not  definitive  nor 
quantitative.  However,  these  data  are  good  indicators  of  the  ability  to 
compensate  for  the  electron  be2ur.  loading,  and  this  is  a  valid  technique  to 
modulate  the  electron  beam. 

As  stated  above,  first  a  cold  test  calibration  was  completed,  determining  the 
cold  resonant  frequency,  the  transfer  function  for  the  cavity  E-dot  probe 
against  the  waveguide  loop  coupler  monitoring  the  magnetron,  and  the 
relationship  between  the  E-dot  probe  and  the  modulation  gap  voltage. 

Then  a  series  of  electron  beam  pulses  where  the  magnetron  frequency  was  varied 
by  1-  or  2-  MHz  steps  from  the  cold  cavity  frequency.  Observations  included 
the  following: 

1)  A  correlation  between  the  time  of  the  reduction  of  the  reflected 
magnetron  power  during  the  time  of  the  electron  beam  pulse  (Fig.  22). 

2)  An  observation  of  signals  from  a  line  of  the  B-dot  array.  The  signal 
amplitude  varied  with  axial  position.  The  results  are  not  quantitative 
because  the  probe  areas  were  not  known  at  the  time  (Fig.  23) .  The  data  do 
show  that  the  double  stub  tuner  allows  coupling  of  rf  power  into  the  cavity. 

This  series  of  experiments  indicates  a  viable  technique  to  modulate  the 
electron  beam  and  that  the  electron  beam  can  carry  the  modulation  to  the  idler 
cavity . 
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Figure  22  Overcoupled  data  showing  reduction  in  reflected  power  during  the  duration  of  the  electron 
beam  pulse. 
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5.0  FUTURE  PIANS 


There  are  three  main  projects  which  will  be  conducted  in  the  near  future: 
validation  of  the  diagnostics;  adding  the  idler  cavity  to  validate  the 
calculated  beam  modulation  coefficient  and  beam  propagation;  and  the  design 
and  construction  of  an  efficient  rf  extractor. 

5.1  VALIDATION  OP  DIAGNOSTICS 

The  various  E-dot  and  B-dot  probes  have  been  independently  calibrated, 
however,  the  interpretation  of  the  data  is  yet  to  be  validated.  For  instance 
in  Section  4  we  discussed  the  radial  E-dot  probe  data  as  an  indicator  for  the 
modulating  gap  voltage.  However,  our  only  cross  check  is  to  measure  the  beam 
modulation  as  the  beam  propagates  past  the  modulating  gap.  Later  on  it  may  be 
useful  to  try  and  use  a  streak  ceimera  to  validate  the  depth  of  the  beeun 
modulation  at  various  axial  positions. 

5.2  IDLER  CAVITY 

The  addition  of  the  idler  cavity  requires  an  accurate  determination  of  the 
axial  drift  distance  to  reach  the  peak  of  the  beam  modulation.  This  axial 
position  may  allow  a  flexibility  of  1  or  2  cm  in  the  exact  position  of  the 
idler  gap.  The  other  question  to  be  addressed  is  the  formation  of  a  virtual 
cathode  in  the  idler  gap  due  to  the  larger  currents  present  in  the  modulated 
electron  beam,  and  the  larger  vacuum  wall  radius.  There  are  the  options  of 
adjusting  the  extent  of  the  axial  gap  and  tuning  the  idler  cavity  by  the 
addition  of  volume  displacements.  These  adjustments  must  be  done  self- 
consistently  to  achieve  the  largest  beam  modulation  at  the  extractor  position. 


5.3  THE  Rf  EXTRACTOR 

The  rf  extractor  technique  is  presently  undefined.  Only  certain  axial 
positions  allow  one  to  extract  the  maximum  rf  power,  due  to  the  phase  bunching 
of  the  electron  beam.  The  research  plan  is  to  begin  with  the  design  of  Fazio 
(Fig.  24) .  Again  the  idea  is  to  validate  or  improve  calculational  techniques 
that  one  might  employ  in  other  HPM  sources.  We  will  simulate  this  rf 
extractor  with  several  different  computational  tools,  such  as  MAGIC  and  finite 
element  codes  like  Hewlett-Packard's  High  Frequency  Structure  Simulator  code. 
Each  type  of  code  brings  different  physics  to  the  problem  and  allows  one  to 
focus  on  special  issues.  We  will  be  looking  at  the  modulated  beeim  electron 
distribution  function,  the  voltage  gen.rrat-J  tl...  if  w.ctractor,  the  power 
coupled  out  by  the  extractor,  and  the  difference  in  the  input  electron  beam 
power  and  the  remaining  beam  power  following  the  rf  extractor.  Again,  one  can 
see  that  self-consistency  checks  are  used  whenever  possible  to  ensure  the 
accuracy  of  the  calculations.  These  are  the  same  checks  that  will  be  used 
experimentally  if  possible.  The  problem  to  be  addressed  has  best  been  stated 
by  Dr.  Keith  Kato:  "Anyone  can  modulate  an  electron  beam,  but  no  one  has  been 
able  to  extract  rf  power."  Finally,  the  rf  extractor  must  incorporate  an 
electron  beam  dump  which  does  not  cause  problems  in  the  rf  extractor 
efficiency  such  as  secondary  electrons  or  plasma  formation. 

This  research  project  has  identified  several  limitations  in  the  way  HPM 
sources,  specifically  amplifiers,  are  developed.  These  limitations  are  being 
addressed  for  future  HPM  source  requirements  so  that  one  may  more  quickly  move 
from  a  concept  or  requirement  to  an  operational  source.  One  limitation  that 
deserves  additional  comment  is  the  material  properties  of  the  metal  used  to 
construct  the  microwave  circuit.  Plasma  production,  secondary  electron 
emission,  multipactoring,  etc.  are  unavailable  qualities  of  the  conductors 
used  in  the  simulations.  Possibly  these  qualities  will  be  included  in  future 
codes  which  will  use  massively  parallel  computers. 
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APPENDICES 


A.  Detsrmlnation  of  B-dot  prob«  aroas 

The  calibration  technique  employed  to  determine  the  effective  area  of  any 
probe  must  include  consideration  of  the  nature  of  the  result  desired.  For 
example,  when  calibrating  B-dot  probes  one  may  wish  information  as  either  a 
current  diagnostic  or  microwave  power  diagnostic.  At  first  glance  one  might 
believe  that  these  calibrations  would  yield  the  same  result,  but  it  will  be 
shown  that  the  nature  of  the  data  obtained  makes  the  calibration  techniques 
differ.  Please  note  that  the  following  calculations  use  MKS  units. 


First,  consider  a  B-dot  probe  used  to  couple  a  microwave  signal  to  a  crystal 
diode  detector.  The  crystal  detector  provides  the  time  variation  of  the  root 
mean  square  (RMS)  power  collected  by  the  probe.  For  the  time  being  ignore  the 
fact  that  typically  some  amount  of  attenuation  is  required  to  keep  from 
overdriving  the  crystal.  The  RMS  power  coupled  to  a  coaxial  cable  is: 


'  ua 


2  Z, 


USE 


(1) 


where  V  is  the  voltage  induced  by  the  B-dot  loop  and  Z^ine  is  typically  50Q. 
The  factor  of  H  (which  is  the  RMS  time  average  factor)  is  the  crucial 
difference  in  the  calibration  analysis. 


The  V  is  related  to  the  Bg  by  Faraday's  Law: 


^  dt 


(2) 


where  Ap„b*  is  the  parameter  to  be  found,  and  Bg  is  the  magnetic  field 
generated  by  the  "Tin  Man"  calibration  line  when  powered  by  continuous  wave 
(CW)  oscillator  at  a  frequency  f=a)/(2  ;i)  . 


One  may  relate  Bg  to  the  current  (I)  or  power  (Pin)  used  to  drive  "Tin  Man"  as: 
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where  Z^in  man  also  50£2,  b  is  radius  of  the  outer  conductor. 
Now  solving  the  above  equations  for  Apnobe* 
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where  Pout/Pin  be  easily  determined  with  a  network  analyzer  such  as  an 
Hewlett-Packard  8510.  The  data  are  obtained  as  attenuation  versus  frequency, 
and  the  area  may  be  determined  from  Equation  5. 

Equation  2  also  allows  determination  of  the  coefficient  between  the  voltage 
generated  on  a  50-0  coaxial  cable  by  a  B-dot  probe  due  to  passage  of  a  current 
pulse.  This  current  pulse  may  be  generated  by  either  an  electron  beaun  or  a 
•short  shot*  from  a  pulse  source.  Essentially  this  technique  is  looking  for 
the  time  variation  of  a  current  pulse,  and  employs  observation  of  the  time 
derivative  of  the  current  pulse  as  follows: 

—  =  c  y(f)  (6) 

dt 


where  C  is  a  coefficient  which  depends  on  the  probe  area  and  position.  Using 
Ampere's  Law  with  Faraday's  Law  one  finds  that  from  C  one  may  get  Aprob.  as 
follows : 


where  C  is  determined  by  comparison  with  other  current  pulse  diagnostics  (i.e. 
current  viewing  resistors  or  previously  calibrated  Rogowskii  coils)  and  for 
this  application  r  =  b  as  in  Equation  5. 

Now  one  would  like  to  find  that  Equation  5  and  Equation  7  yield  the  same 
answers  for  the  effective  area  for  the  B-dot  probes.  However,  the  ratio  of 
Equation  7  to  Equation  5  is  found  to  equal  V2.  This  difference  is  the  result 
of  the  nature  of  the  data  provided  by  the  B-dot  probe. 


B.  SUPERPISH  INPUT 


Data  file  for  the  RKA  geometry,  input  to  AUTOMESH: 

X  RKA  SINGLE  CAVITY (3/24/92) 

$REG  NREG=4,  DX=0.5,  DY=0.5.  XMAX=45.0.  yMAX=15.0,  MAT=1,  NP0INT=5, 

NDRIVE=1,  IB0UND=1,  $END 
$P0  X=0.0,  Y=0.0,  $END 
$P0  X=0.0,  Y=15.0,  $END 
$P0  X=45.0,  Y=15.0,  $END 
$P0  X=45.0,  Y=0.0,  $END 
$P0  X=0.0,  Y=0.0,  $END 
$REG  NP0INT=11,  MAT=0,  $END 
$P0  X=0.0,  Y=7.5,  $END 
$P0  X=12.0,  Y=7.5,  $END 
$P0  X=12.0,  Y=14.5,  $END 
$P0  X=25.0,  Y=14.5,  $END 
$P0  X=25.0,  Y=8.5,  $END 
$P0  X=13.5,  Y=8.5,  $END 
$P0  X=13.5,  Y=7.5,  $END 
$P0  X=45.0,  Y=7.5,  $END 
$P0  X=45.0,  Y=15.0,  $END 
$P0  X=0.0,  Y=15.0,  $END 
$P0  X=0.0,  Y=7.5,  SEND 
$REG  NP0INT=8,  MAT=0,  SEND 
SPO  X=0.0,  Y=0.0,  SEND 
SPO  X=0.0,  Y=1.0,  SEND 
SPO  X=2.0,  Y=1.0,  SEND 
SPO  X=6.0,  Y=6.5,  SEND 
SPO  X=:6.0,  Y=6.0,  SEND 
SPO  X=2.5,  Y=1.0,  SEND 
SPO  X=2.5,  Y=0.0,  SEND 
SPO  X=0.0,  Y=0.0,  SEND 
SREG  NP0INT=1,  SEND 
SPO  X=i2.0,  Y=14.50,  SEND 

CON  input  for  LATTICE: 

*19  1  *21  1  0  0  0  s 

This  sets  cylindrical  symmetry,  and  allows  for  Neumann  boundary  conditions  on 
the  left  and  right  surfaces  of  the  problem.  This  choice  is  required  when  no 
conducting  surface  matches  the  left/right  edges. 

PSFPLOT  input  to  plot  geometry  and  mesh: 

in  response  to  (num,  itri,  nphi,  inap,  nswxy) 

0,0, 0,0,0  to  plot  problem  geometry,  and 

0,1, 0,0,0  to  plot  triangular  mesh  used  to  solve  problem. 


CON  input  for  : 

*65  initial  frequency  guess  to  find  a  specific  mode  s; 
*62  nsteps  *63  Ak^  *65  initial  f requency ( fo)  s. 
fj  -  final  frequency  of  search  range 
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At* 


c  (nsteps  -  1) 


PSFPLOT  input  to  plot  electric  field  contours: 
in  response  to  (num,  itri,  nphi,  inap,  nswxy) 

1,0, #,0,0  where  #  specifies  the  number  of  electric  field  contours  to  plot. 
Input  for  SHY  to  specify  normalized  electric  field  amplitude: 

*42  kmin  *43  ktop  *44  Imin  *45  Imax  *74  ascale 

where  kmin/kmax  and  Imin/lmax  specify  the  left/right  axial  and  lower/upper 
radial  mesh  points  for  the  electric  field  amplitude.  Notice  that  if  lmin=lmax 
then  the  normalized  Ej  is  found  along  the  contour  from  kmin  to  kmax.  Also,  if 
kmin=kmax  then  the  normalized  E^  is  found  along  the  contour  from  Imin  to  Imax. 
One  must  then  select  the  proper  value  from  the  listed  normalized  field 
amplitudes  for  the  desired  mesh  coordinate.  The  parameter  ascale  may  be  set 
to  1.0  if  the  absolute  value  is  not  required,  or  any  other  normalizing  value 
may  be  used  to  properly  scale  the  electric  field. 


39 


C.  MAGNET  INPUT 


MAGNET  was  originally  written  by  John  Freeman  of  Sandia  National  Laboratory. 
The  code  was  later  modified  by  several  users  as  the  code  was  ported  to  other 
platforms  such  as  PCs. 

One  may  use  CDFMAG  to  input  the  geometry  and  current  for  each  winding  in  the 
coil,  or  write  a  short  FOPTRAN  program  to  generate  the  input  deck.  The  format 
of  each  line  is  critical,  but  is  documented  within  the  code.  Basically,  the 
format  is: 

#;  for  the  total  number  of  lines  in  the  data  file 

r(in.),  left  z(in.),  right  2(in.),  #  turns,  I  (Amps);  this  line  is  repeated  # 
times.  The  code  MAGNET  does  require  the  inputs  of  dimensions  in  inches,  this 
is  a  minor  irritant  and  has  not  been  updated.  The  following  geometry  list  is 
from  a  short  Fortran  program.  Other  techniques  are  possible  to  create  the 
file,  as  long  as  the  format  of  the  line  is  followed. 

For  the  RKA  magnets  #=550 

Sample  MAGNET  input  for  the  RKA  coils: 

550 

Coil  1 

8.737500E+00  6.560000E+00  9.B60000E+00  1  l.OOOOOOE+00 

8.762500E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.787499E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.812499E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.837499E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.862498E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.887498E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.912498E+00  6.560000E+00  9.560000EtOO  1  l.OOOOOOE+00 

8.937497E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

8.962497E+00  6.560000E+00  9 . 560000E+00  1  l.OOOOOOE+00 

8.987496E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.012496E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.037496E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.062495E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.087495E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.112494E+00  6.560000E+00  9.560000Et00  1  l.OOOOOOE+00 

9.137494E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.162494E+00  6.560000Et00  9.56000CE+00  1  l.OOOOOOE+00 

9.187493E+00  6 . 560000E*00  9.560000E+00  1  l.OOOOOOE+00 

9.212493E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.237493E+00  6.560000E+00  9.560000Et00  1  l.OOOOOOE+00 

9.262492E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.287492E+00  6.560000E+00  9.560000Et00  1  l.OOOOOOE+00 

9.312491E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.337491E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.362491E+00  6.B60000E+00  9.S60000EtOO  1  l.OOOOOOE+00 

9.387490E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.412490E+00  6.560000E+00  9.560000EtOO  1  l.OCOOOOE+00 

9.437490E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.462489E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.487489E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.512488E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.B37488E*00  6.560000E+00  9.560000E*00  1  l.OOOOOOE+00 

9.562488E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.587487E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.612487E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.637486E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.662486E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.687486E*00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.712485E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9,737485E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.762485E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.787484E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.812484E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.837483E+00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+00 

9.862483Et00  6.560000E+00  9 . 560000E+ 00  1  l.OOOOOOE+00 

9.887483Et00  6.560000E+00  9.560000E+00  1  l.OOOOOOE+OO 


9.912482E+00 
9.937482E+00 
9.962481E+00 
9.987481E+00 
1.001248E+01 
1.003748E+01 
1 .006248E+01 
1.008748E+01 
1.011248E+01 
1.013748E+01 
I . ji6^4S£t01 
1.018748E+01 
1.021248E+01 
1.023748E+01 
1.026248E+01 
1.028748E+01 
1.031248E+01 
1.033748E+01 
1.036248E+01 
1.038748E+01 
1.041247E+01 
1.043747E+01 
1.046247E+01 
1.048747E+01 
1.051247E+01 
1.053747E+01 
1.056247E+01 
1.058747E+01 
1.061247E+01 
1 .063747E+01 
1.066247E+01 
1.068747E+01 
1.071247E+01 
1.073747E+01 
1.076247E+01 
1.078747E+01 
1.081247E+01 
1.083747E+01 
1.086247E+01 
1.088747E+01 
1.091247E+01 
1.093747E+01 
1.096247E+01 
1.098747E+01 
1.101247E+01 
1.103747E+01 
1.106246E+01 
1.108746E+01 
1.111246E+01 
1.113746E+01 
1.116246Et01 
1 . 118746E+01 
1.121246E+01 
1.123746E+01 
1.126246E+01 
1 . 128746E+01 
1.131246E+01 
1.133746E+01 
1.136246E+01 
1.138746E+01 
1 .141246E+01 
1.143746E+01 
1.146246E+01 
1 . 148746E+01 
1.151246E+01 
1.153746E+01 
1.156246E+01 
1.158746E+01 
1.161246E+01 
1.163746E+01 
1 . 166246E+01 
1 . 168746E+01 
1.171245E4.01 
1 . 173745E+01 
1 . 176245Et01 
1 . 178745E+01 
1 . 181245E4.01 
1 . 183745Et01 


6.560000E+00 
6.560000E+00 
6.560000E+00 
6.560000E+00 
6 .560000E+00 
6.560000E+00 
6 .560000E+00 
6.560000E+00 
6 .560000E+00 
6.560000E+00 
£ .5607' OE+00 
6.560000E+00 
6 .560000E+00 
6.560000E+00 
6.560000E+00 
6 .560000E+00 
6 .560000E+00 
6.B60000E+00 
6 .560000E+00 
6 .560000E+00 
6 .560000E+00 
6.560000E+00 
6.560000E+00 
6.560000E+00 
6 .560000E+00 
6.560000E+00 
6.560000E+00 
6 .560000E+00 
6.560000E+00 
6 .560000E+00 
6.B60000E+00 
6.B60000E+00 
6 .B60000E+00 
6.560000E+00 
6.B60000E+00 
6.B60000E+00 
6.560000E+00 
6.560000E+00 
6 .560000E+00 
6.560000E*00 
6.560000E+00 
6.560000E+00 
6.B60000E+00 
6  .B60000E+00 
6.B60000E+00 
6.560000E+00 
6.B60000E+00 
6.B60000E+00 
6.B60000E+00 
6.560000E+00 
6.B60000E+00 
6.B60000E+00 
6.B60000E+00 
6.560000E+00 
6.560000E+00 
6  .B60000E+00 
6  .B60000E+00 
6  .B60000E+00 
6 .B60000E+00 
6  .B60000E+00 
6.B60000E+00 
6 .560000E+00 
6.B60000E+00 
6 .B60000E+00 
6  .B60000E+00 
6.B60000E+00 
6.B60000E+00 
6.B60000E+00 
6.560000E+00 
6 .560000E+00 
6.B60000E+00 
6.B60000E+00 
6 .B60000E+00 
6 .B60000E+00 
6 .B60000E+00 
6.B60000E+00 
6 . 560000E+00 
6 .560000E*-00 


9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.560000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.560000E+00 
9.B60000E+00 
9.B60000E+00 
7 .  *^o0r»r»r»E-^ GO 
9.560000E+00 
9 .B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.560000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E400 
9.560000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E+00 
9.B60000E400 
9.B60000E400 
9.B60000E400 
9.B60000E400 
9.S60000E400 
9.560000E400 
9.560000E400 
9.B60000E400 
9.560000E400 
9.B60000E400 
9.S60000E400 
9.S60000E400 
9.560000E400 
9.B60000E400 
9.560000E400 
9.560000E400 
9.S60000E400 
9.B60000E400 
9.S60000E400 
9.B60000E400 
9 .560000E400 
9.560000E400 
9.B60000E400 
9  .B60000E400 
9.560000E400 
9.B60000E400 
9.560000E400 
9.560000E400 
9.560000E400 
9.B60000E400 
9.560000E400 
9.560000E400 
9.B60000E400 
9.560000E400 
9.560000E400 
9.560000E400 
9.560000E400 
9.560000E400 
9.B60000E400 
9.560000E400 
9 .560000E400 
9.560000E400 
9.560000E400 
9.560000E400 
9.B60000E400 
9 .560000E400 
9.560000E400 
9.B6OOOOE4OO 
9 .560000E400 
9.560000E400 


1  I.OOOOOOE4OO 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  l.OOOOOOE+00 
1  1.000000E400 
1  1.000000E400 
*  I.GGCGOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 

1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  l.OOOOOOE+00 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
I  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 

1  1.000000E400 

1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  1.000000E400 
1  I.OOOOOOE4OO 
1  I.OOOOOOE4OO 
1  1.000000E400 
1  1.000000E400 
\  1.000000E400 
1  I.OOOOOOE4OO 
1  1.000000E400 
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1.18624BE+01  6.560000E+00 

1.188745E+01  6.560000E+00 

1.191245E+01  6.560000E+00 

1.193745E+01  6.560000E+00 

1.196245E+01  6.560000E+00 

1.198745E+01  6.B60000E+00 

1.20124BE+01  6.B60000E+00 

1.20374BE+01  6.B60000E+00 

1.20624BE+01  6.B60000E+00 

1.20874BE+01  6.B60000E+00 

:i  e .B6''o"9r+oo 

1.21374BE+01  6.B60000E+00 

1.216245E+01  6.B60000E+00 

1.21874BE+01  6.B60000E+00 

1.22124BE+01  6.B60000E+00 

1.22374BE+01  6.B60000E+00 

1.22624BE*01  6.B60000E+00 

1.22874BE+01  6.B60000E+00 

1.23124BE+01  6.B60000E+00 

1.23374SE+01  6.S60000E+00 

1.236244E+01  6.B60000E+00 

1.238744E+01  6.B60000E+00 

1.241244E+01  6.B60000E+00 

1.243744E+01  6.B60000E+00 

1.246244E+01  6.B60000E+00 

1.248744E+01  6  .  B60000E*-00 

1.2B1244E+01  6.B6OOO0EtOO 
1.2B3744E+01  6 . B60000E+00 

1.2B6244E+01  6.B60000E+00 

1.2B8744E+01  6.B60000E»00 

1.261244E+01  6.B60000E+00 

1.263744E+01  6.B60000EtOO 

1.266244E+01  6.B60000E+00 

1.268744E+01  6.B60000E+00 

1.271244E+01  6.B60000E+00 

1.273744E+01  6.B60000E+00 

1.276244E+01  6.B60000E+00 

1.278744Et01  6.B60000E+00 

1.281244E+01  6.B60000E+00 

1.283744E+01  6.B60000E+00 

1.286244E+01  6.B60000E*00 

1.288744E+01  6.S60000EtOO 

1.291244E+01  6.B60000E+00 

1.293744Et01  6.B60000E+00 

1.296244E+01  6.S60000E+00 

1.298744E+01  6.S60000E+00 

1.301243E+01  6.B60000E+00 

1.303743E+01  6 . B60000E+00 

1.306243E+01  6.S60000E+00 

1.308743E+01  6.560000E*00 

1.311243E+01  6.B60000Et00 

1.313743E+01  6  .■56OOOOE+OO 

1.316243E+01  6,B60000E+00 

1.318743E+01  6  B60000E+00 

1.321243E+01  6.560000E»00 

Coil  2 

8.737B00Et00  1.B36000E+01 

8.762S00E+00  1.B36000E+01 

8.787499E+00  1.B36000E+01 

8.812499E+00  1.536000E+01 

8.837499E+00  1.B36000E+01 

8.862498E+00  1.B36000E+01 

8.887498E+00  1,S36000E+01 

8.912498E+00  1.B36000E+01 

8.937497E+00  1.B36000E+01 

8 . 962497E't00  1.S36000E+01 

8.987496E+00  1.B36000E+01 

9.012496E+00  1.B36000E+01 

9.037496E+00  1 .  B36000E-^01 

9.06249BE+00  l.B36000Et01 

9.08749BE+00  1.B36000E+01 

9.112494E+00  l.B36000E*01 

9.137494E+00  1.B36000E+01 

9.162494E+00  1.B36000E+01 

9.187493E+00  1.536000E+01 

9.212493E+00  l.B36000E*01 

9.237493E*00  1  .  B3 6000E+ 0 1 


9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E-t00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+CO  1  1 .OOOCOOEi 00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  1 . OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.S60000E<00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.S60000E+00  1  l.OOOOOOE+00 
9.B60000Et00  1  l.OOOOOOE+00 
9.B60000Et00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E>00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.S60OOOE-fOO  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000Et00  1  1.000000E*00 
9.S60000Et00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E*00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.660000E<-00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000Et00  1  l.OOOOOOE+00 
9.B60000Et00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 
9.B60000E+00  1  l.OOOOOOE+00 


1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000Et01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  i.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000Et01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836n00E*01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOEtOO 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000EtCl  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
1.836000E+01  1  l.OOOOOOE+00 
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9.262492E+00 
9.287492E+00 
9.312491E+00 
9.337491E+00 
9.362491E+00 
9.387490E+00 
9.412490E+00 
9.437490E+00 
9.462489E+00 
9.487489E+00 
9.512488E+00 
9.5374bafc+00 
9.562488E+00 
9 .587487E+00 
9.612487E+00 
9.637486E+00 
9.662486E+00 
9.687486E+00 
9.712485E+00 
9.737485E+00 
9.762485E+00 
9.'’87484E+00 
9.812484E+00 
9.837483E+00 
9.862483E+00 
9.887483E+00 
9.912482E+00 
9.937482E+00 
9.962481E+00 
9.987481E+00 
1.001248E+01 
1 .003748E+01 
1.006248E+01 
1.008748E+01 
1.011248E+01 
1 .013748E+01 
1.016248E+01 
1.018748E+01 
1.021248E+01 
1.023748E+01 
1.026248E+01 
1.028748E+01 
1 .031248E+01 
1.033748E+01 
1.036248E+01 
1.038748E+01 
1 .041247E+01 
1.043747E+01 
1.046247E+01 
1.048747E+01 
1.051247E+01 
1 .053747E+01 
1 .056247E+01 
1.058747E+01 
1.061247E+01 
1 .063747E+01 
1.066247E+01 
1.068747E*01 
1.071247E+01 
1 .073747E+01 
1.076247E+01 
1.078747E+01 
1.081247E+01 
1.083747E+01 
1 .086247E+01 
1.088747E+01 
1 .091247E+01 
1 .093747E+01 
1 .096247E+01 
1.098747E+01 
1 . 101247E+01 
1.103747E+01 
1.106246E+01 
1.108746E+01 

Coil  3 

8 .737500E+00 
8 .762500E*00 


1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
i.5360v>0E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1.536000E+01 
1.S36000E+01 
1 .536000E+01 
1.536000E+01 
1 .536000E+01 
1 .536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.B36000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1 .536000E+01 
1.536000E+01 
1 .536000E+01 
1.B36000E+01 
1.536000E+01 
1.536000E+01 


2 .286000E+01 
2 .286000E+01 


1.836000E+01  i 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1 . 8JbD00i:,+  01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E<-01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836700E+01  1 
1.836000E+01  1 
1,836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000Et01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000Et01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000Et01  1 
1 .836000E*01  1 
1.836000E+01  1 
1.836000E+01  1 
1.836000E+01  1 


2 .586000Et01  1 

2.586000E+01  1 


l.OOOOOOE+00 
l.OOOOOOE+00 
1 . ooooooE+o: 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
1 .OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
1  -OOOoOCjb+uO 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
X.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 
1 .OOOOOOE+00 
l.OOOOOOE+00 
1 .OOOOOOE+00 


1 .OOOOOOE+00 
1 .OOOOOOE+00 
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8 .787499E+00 
8.812499E+00 
8.837499E+00 
8.862498E+00 
8.887498E+00 
8 . 912498E+00 
8.937497E+00 
8.962497E+00 
8.987496E+00 
9.012496E+00 
9.037496E+00 
'i .  'Joi4abbtOu 
9.087495E+00 
9 . 112494E+00 
9.137494E+00 
9.1G2494E+00 
9.187493E+00 
9.212493E+00 
9.237493E+00 
9.262492E+00 
9.287492E+00 
9.312491E+00 
9.337491E+00 
9.362491E+00 
9.387490E+00 
9.412490E+00 
9.437490E+00 
9.462489E+00 
9.487489E+00 
9.512488E+00 
9.537488E+00 
9.B62488E+00 
9.587487E+00 
9.612487E+00 
9.637486E+00 
9.662486E+00 
9.687486E+00 
9.712485E+00 
9.737485E+00 
9 .762485E+00 
9.787484E+00 
9 .812484E+00 
9.837483E+00 
9.862483E+00 
9.887483E+00 
9.912482E+00 
9.937482E+00 
9.962481E+00 
9.987481E+00 
1.001248E+01 
1.003748E+01 
1.006248E+01 
1.008748E+01 
1.011248E+01 
1.013748E+01 
1.016248E+01 
1.018748E+01 
1.021248E+01 
1 .023748E+01 
1 .026248E+01 
1 .028748E+01 
1 .031248E+01 
1.033748E+01 
1 .036248E+01 
1.038748E+01 
1 .041247E+01 
1.043747E+01 
1 .046247E+01 
1 .048747E+01 
1 .051247E+01 
1.053747E+01 
1.056247E+01 
1 .058747E+01 
1.06i247E+01 
1 .063747E+01 
1.066247E+01 
1 .068747E+01 
1.071247E+01 


2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 

2  .ZOOUUUCj.'V.  X 

2.286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
;  .286000E+01 
2.286C30E+01 
2.2860C0E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .2e6000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000Et01 
2.2OG000E+01 
2 .286000E+Oi 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2.286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2  .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 
2 .286000E+01 


2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 

^  ATT*  ^  <4  ^ 

2 .586000Et01 
2.586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2.586000E+01 
2 .586040E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2.586000E+01 
2.586000E+01 
2 .586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
■  5»«0n0E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E*-01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2.586000E+01 
2 .B86000E+01 
2 .586000E+01 
2 .586000E+01 
2.586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E*01 
2 .586000E*01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 . 586000E+01 
2 .586000E+01 
2.586000E401 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E+01 
2 .586000E401 
2.586000E+01 


1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  1  oooooriE*on 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+00 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 

1  i.oooooop+c: 

1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOEiOO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOCE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+CO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
1  l.OOOOOOE+OO 
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1.073747E+01  2.286000E+01 
1.076247E+01  2.286000E+01 
1.078747E+01  2.286000E+01 
1.081247E+01  2.286000E+01 
1.083747E+01  2.286000E+01 
1.086247E+01  2.286000E+01 
1.088747E+01  2.28d000E+01 
1.091247E+01  2.286000E+01 
1.093747E+01  2.286000E+01 
1.096247E+01  2.286000E+01 
1.098747E+01  2.286000E+01 
1.101247E+01  2.286000E+01 
1.103747E+01  2.286000E+01 
1.106246E+01  2.286000E+01 
1.108746E+C1  2.286000E+01 

Coil  4 

8.737500E+00  3.156000E+01 
8.762500E+00  3.156000E+01 
8.787499E+00  3.156000E+01 
8.812499E+00  3.156000E+01 
8.837499E+00  3.156000E+01 
8.862498E+00  3.156000E+01 
8.887498E+00  3.156000E+01 
8.912498E+00  3.156000E+01 
8.937497E+00  3.156000E+01 
8.962497E+00  3.156000E+01 
8.967496E+00  3.156000E+01 
9.012496E+00  3.156000E+01 
9.037496E+00  3.156000E+01 
9.062495E+00  3.156000E+01 
9.087495E+00  3.156000E+01 
9.112494E+00  3.156000E+01 
9.137494E+00  3.156000E+01 
9.162494E+00  3.156000E+01 
9,187493E+00  3.156000E+01 
''.212493E+00  3.156000Et01 
9.23?4?3F+nO  3.156000E+01 
9.262492E+00  3  156000E401 
9.2R7492E+00  3  . 156u0i/E+01 
9.312491E+00  3.156000E+01 
9.337491E+00  3.156000E+01 
9.362491E+00  3.1S6000E+01 
9.387490E+00  3.156000E+01 
9.412490E+00  3.156000E+01 
9.437490E+00  3.156000E+01 
9.462489Et00  3.156000E+01 
9.487489E+00  3.156000E+01 
9.B12488E+00  3.156000E+01 
9.B37488E+00  3.156000E+01 
9.562488E+00  3.156000E+01 
9.587487E+00  3.156000E401 
9.612487E+00  3.156000E+01 
9.637486E+00  3.156000E+01 
9.662486Et00  3.156000E+01 
9.687486E+00  3.156000E+01 
9.712485E+00  3.156000E+01 
9.737485E+00  3.156000E+01 
9.762485E+00  3.156000E+01 
9.787484E+00  3.156000E+01 
9.812484E+00  3.156000E+01 
9.837483E+00  3.156000E+01 
9.862483E+00  3.156000E+01 
9.887483E+00  3.1B6000E+01 
9.912482E+00  3.156000E+01 
9.937482E+00  3.156000E+01 
9.962481E+00  3.156000E+01 
9.987481E+00  3.156000E+01 
1.001248E+01  3.156000E+01 
1.003748E+01  3.156000E+01 
1.006248E+01  3.156000E+01 
1.008748E+01  3.156000E+01 
1.011248E+01  3.156000E+01 
1.013748E401  3.156000E+01 
1.016248E401  3.156000E401 
1.018748E+01  3.156000E401 
1.021248E+01  3.156000E+01 
1.023748E401  3.156000E+01 


2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.B86000E+01  1  l.OOOOOOE+00 
2.586000E+01  1  l.OOOOOOE+00 
2.B86000E401  1  l.OOOOOOE+00 
2.B86000E+01  1  l.OOOOOOE+00 
2.B86000Et01  1  l.OOOOOOE+00 
2.B86000E+01  1  l.OOOOOOE+00 


3.4B6000Et01  1  l.OOOOOOE+00 
3.4B6000E401  1  l.OOOOOOE+00 
3 .4B6000E401  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  i  l.OOOOOOE+00 
3  4S6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4S6000Et01  1  l.OOOOOOE+00 
3.4B6000E+01  1  1  OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E401  I  l.OOOOOOE+00 
3.4B6000Et01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4S6U00E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.4B6000E+0i  1  l.OOOOOOE+00 
3.4B6000Et01  1  l.OOOOOOE+OO 
3.4B6000E401  1  l.OOOOOOE+OU 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4S6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3 .4‘^6000E+01  1  i.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000Et01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
3.4B6000E+01  1  l.OOOOOOE+OO 
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1.026248E*01  3.156000E+01 
1.028748E+01  3.156000Et01 
1.031248E+01  3.156000E+01 
1.033748E+01  3.156000E+01 
1.036248E+01  3.156000E+01 
1.038748E+01  3.156000E+01 
1.041247E+01  3.156000E+01 
1.043747E+01  3.156000E+01 
1.046247E+01  3.156000E+01 
1.048747E+01  3.156000E+01 
1.051247E+01  3.156000E+01 
i  .  03  3  74  7^1- ox  C.ljoCOOExOl 
1.056247Et01  3.156000E+01 
1.058747E+01  3.156000E+01 
1.061247E+01  3.156000E+01 
1.063747E+01  3.156000E+01 
1.066247E+01  3.156000Ex01 
1.068747E+01  3.156000E+01 
1.071247E+01  3.156000E+01 
1.073747E+01  3.156000E+01 
1.076247E+01  3.156000E+01 
1.078747Ex01  3.156000E+01 
1.081247E+01  3.156000E+01 
1.083747E+01  3.156000E+01 
1.086247E+01  3.1B6000E+01 
1.088747E+01  3.156000E+01 
1.091247E+01  3.156000E+01 
1.093747E+01  3.156000E+01 
1.096247E+01  3.1S6000E+01 
1.098747E+01  3.15COOOEt01 
1.101247E+01  3.156000E4-01 
1.103747E+01  3.156000E+01 
1.106246E+01  3.156000E+01 
1.108746E+01  3.156000Ex01 
1.111246E+01  3.156000Ex01 
1.113746E+01  3.156000Ex01 
1.116246E+01  3.:56000E+01 
1.118746E+01  3.156000E+01 
1.121246E+01  3.156000EX01 
1.123746E+01  3.:56000E+01 
1.126246E+01  3.156000E+01 
1.128746Ex01  3.156000E+01 
1.131246E+01  3.156000EX01 
1.133746E+01  3.156000E+01 
1.136246E+01  3.156000E+01 
1.138746E+01  3.156000Ex01 
l.i4:246E+01  ?  156000Et01 
1.143746E+01  3.156000E+01 
1.146246E+01  3.i56000Ex01 
1.14e746E+01  3.156000E+01 
1.151246E+01  3.156000E+01 
1.153746E+01  3.156000E*01 
1.156246E+01  3.156000E+01 
1.158746Ex01  3.156000Ex01 
1.161246E+01  3.156000E+01 
1.163746E+01  3.156000Et01 
1 .  :‘;6246E+01  3.156000E+01 
1.168746E+01  3.156000E+01 
1.171245E+01  3.156000Ex01 
1.173745E+01  3.156000E+01 
1.176245E+01  3.156000E*01 
1.178745E+01  3.156000E+01 
1.181245E+01  3.156000E+01 
1.183745Et01  3.1B6000E+01 
1.186245E+01  3.156000E+01 
1.188745E+01  3.156000E+01 
1.191245Ex01  3.156000E+01 
1.193745E+01  3.1B6000E+01 
1.196245Ex01  3.156000Ex01 
1.198745E+01  3.156000E+01 
1.201245EX01  3.1B6000E+01 
1.203745E+01  3.156000E+01 
1.206245Ex01  3.156000E+01 
1.208745Ex01  3.156000E+01 
1.211245E+01  3  ,lB6000E-f01 
1.213745E+01  3.156000Ex01 
1.216245E+01  3.156000E+01 
1.218745E+01  3.156000E+01 


3,456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l,000000E+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000ExQl  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000EX01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000Ex01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E*01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOExOO 
3.456000E+01  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  1 . OOOOOOEt 00 
3.456000E+01  1  l.OOOOOOE+00 
3.456000E401  1  l.OOOOOOE+00 
3.456000E+Oi  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOEtOO 
3.456000E401  1  l.OOOOOOE+OO 
3.456000E401  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+C1  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOCE+00 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.4B6000Et01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3 .456C''OE+01  1  l.OOOOOOE+OO 
3.456000E+C1  1  l.OOOOOOE+00 
3.4B6000E+01  1  l.OOOOOOE+00 
3.456000E+01  1  l.OOOOOOE+OO 
3.45600OE+01  1  l.OOCOOOE+00 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
3.456000E+01  1  l.OOOOOOE+OO 
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1 .221245E+01 

3 . 156000E+01 

3 .4B6000Et01 

1 

l.OOOOOOE+00 

1 .223745E+01 

3 . lB6000Et01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1.22624BE+01 

3 . 156000E+01 

3 .4S6000E+01 

1 

l.OOOOOOE+00 

1 .228745E+01 

3.1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1 .231245E+01 

3 . 1B6000E+01 

3  .4B6000E+01 

1 

1  .OOOOOOE+00 

1.233745E+01 

3 . 1B6000E+01 

3  .4S6000E+01 

1 

1  .OOOOOOE+00 

1 .236244E+01 

3  .1B6000E4-01 

3 .4S6000E+01 

1 

l.OOOOOOE+00 

1 .238744E+01 

3.1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1.241244E+01 

3.1B6000E+01 

3 .456000E+01 

1 

l.OOOOOOE+00 

1 .243744E+01 

3  .  1B6000E+01 

3  .4B6000E+01 

1 

1  .OOOOOOE+00 

1.246244E+01 

3.1B6000E+01 

3  .4B6000E+01 

1 

l.OOOOOOE+00 

1.243'’ ■)4Et01 

3 . 1B6000E+01 

3 .4S6000E+01 

1 

1  .OOOOOOE+00 

1 .251244Et01 

3 .1B6000E+01 

3 .4B6000E+01 

1 

1  .OOOOOOE+00 

1.253744E+01 

3 .1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1 .2S6244E+01 

3 . 1S6000E+01 

3 .4S6000E401 

1 

l.OOOOOOE+00 

1 .258744E+01 

3 .1B6000E+01 

3 .4S6000E+01 

1 

1  .OOOOOOE+00 

1 .261244E+01 

3 . 1B6000E+01 

3  .4B6000E*01 

1 

l.OOOOOOE+00 

1.263744E+01 

3  .  1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1.266244E+01 

3.1B6000E+01 

3 .4B6000E+01 

1 

1 .OOOOOOE+00 

1 .268744E+01 

3.1B6000E+01 

3.4B6000E*01 

1 

l.OOOOOOE+00 

1.271244E+01 

3.1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1 .273744E+01 

3  .  1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1.276244E+01 

3.1B6000E+01 

3 .4B6000E+01 

1 

1  .OOOOOOE+00 

1.278744Et01 

3  .1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1 .281244E+01 

3.1B6000E+01 

3  .4B6000E*01 

1 

l.OOOOOOE+00 

1 .283744E+01 

3 . 1B6000E+01 

3 .4B6000Et01 

1 

l.OOOOOOE+00 

1 .286244E+01 

3  .  1B6000E*01 

3 .4B6000E+01 

1 

l.OOOOOOE+00 

1 .288744E+01 

3 . 1B6000E+01 

3 .4S6000Et01 

1 

l.OOOOOOE+00 

1.291244E+01 

3  .  1B6000E+01 

3 .4B6000E+01 

1 

l.OOOOOOE+OC 

1 .293744E+01 

3  .1B6000E+01 

3 .4B6000E+01 

1 

1  .OOOOOOE+00 

1 .296244E+01 

3 . 1B6000E+01 

3  .4B6000E+O1 

1 

1  .OOOOOOE+00 

1.298744E+01 

3  .1B60C0E+01 

3 .4B6000Et01 

1 

l.OOOOOOE+00 

1 .301243E+01 

3.1B6000E+01 

3  .4S6000E+01 

1 

l.OOOOOOE+00 

1 .303743E+01 

3 . 1B6000E+01 

3 .4B5000E+01 

1 

l.OOOOOOE+00 

1.306243E+01 

3  .lB6000Et01 

3  .4S6000E+01 

1 

i.OOOOOOE+00 

1 .308743Et01 

3  .  1B6000E+01 

3 .4S6000E+01 

1 

l.OOOOOOE+00 

1 .311243E*01 

3  .  lB6000Et01 

3  .4B6000E+01 

1 

l.OOOOOOE+00 

1.313743E+01 

3  .  lB6000Et01 

3  .4S6000E+01 

1 

1  .OOOOOOE+00 

1.31S243E+01 

3  .  1B6000E+01 

3  .4S6000E+01 

1 

l.OOOOOOE+00 

1.318743E+01 

3.1B6000E*01 

3  .4B6000E+01 

1 

l.OOOOOOE+00 

1.321.i43E*01 

3 . 1B6000E+01 

3  ,  4S6000E+01 

1 

1 .OOOOOOE+00 

Sample  geometry  input  for  POISSON: 

RKA  MAGNET  (7/1/92) 

$REG  NREG=5 , XMAX=80 . 0 , DX=1 . , YMIN=0 . , YMAX=190 . , NP0INT=5 , $END 

$P0  X=0.0,  Y=0. ,  $END 

$P0  X=80.0,  Y=0.0,  $END 

$P0  X=80.0,  Y=190.0,  SEND 

$P0  X=0.0,  Y=190.0,  SEND 

SPO  X=0.0,  Y=0.0,  SEND 

SREG  MAT=1,  CUR=3.40E5,  NP0INT=5,  SEND 

SPO  X=22.19,  Y=59.10,  SEND 

SPO  X=33.56,  Y=59.10,  SEND 

SPOX=33.56,  Y=66.72,  SEND 

SPO  X=22.19,  Y=66.72,  SEND 

SPO  X=22.19,  Y=59.10,  SEND 

SREG  MAT=1,  CUR=1.80E5,  NP0INT=5,  SEND 

SPO  X=22.19,  Y=81.45,  SEND 

SPO  X=28.16,  Y=81.45,  SEND 

SPO  X=28.16,  Y=89.07,  SEND 

SPO  X=22.19,  Y=89.07,  SEND 

SPO  X=22.19,  Y=81.45,  SEND 

SREG  MAT=1,  CUR=1.80E5,  NPOINT=5,  SEND 

SPO  X=22.19,  Y=100.50,  SEND 

SPO  X=28.16,  Y=100.50,  SEND 

SPO  X=28.16,  Y=108.12,  SEND 

SPO  X=22.19,  Y=108.12,  SEND 
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$P0  X=22.19, 
$REG  MAT=1, 
$PO  X=22. 19, 
$P0  X=33 .56, 
$P0  X=33 .56, 
$P0  X=22.19, 
$P0  X=22.19, 


Y=100.50, 
CUR=3 .40E5 
Y=122.60, 
Y=122.60, 
Y=130 .22, 
Y=130.22, 
Y=122.60, 


$END 

NP0INT=5, 

$END 

$END 

$END 

SEND 

SEND 


SEND 


One  then  goes  through  LATTICE 
except  *19  1  s.  PSFPLOT  is  us 


Sample  2-D  grid  input  for  the 


and  POISSON,  no  extra  CON  variable 
ed  as  in  . 

magnetic  diffusion  code  MAGDIF: 


must  be  set 


Note : 
turns 
coil. 


An  index  of  4  indicates  the  stainless  steel, 
for  the  larger  coil,  an  index  of  6  indicates 
and  an  index  of  0  indicates  free  space. 


an  index  of  5 
the  turns  for 


indicates 
a  smaller 


*aba  outline  3/24/92 
Sinp  nr=120  n2=428  dr=0.0050  d2=0.0050 
npas  =  4  nact=2  sigma=7 . 854 , 1 ., 1 ., 1 . 745, 0 . 2629 , 0 . 3480 ,  1 ., 1 ., 1 . , 
amu=100000. , l.,l.,l.,l.,l.,l.,l.,l. 

Send 


0000000 0000000000 00 00 00 OOOOOOvuUUbOsiyvvvCCOOCOCOCOOCSCOC^OOOOOOOOOO 000000000000 0000000000 00 OOOOOOOOOOOOOOOOOOOOOOOO 00000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOODOOOOOO 

OOOOC 0000000000000000 0000000000000000000000 ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
00000000000000000000000000000 0000 000000000000000000000000000000 000000 000000 ooooooooooooooooooooooooooooooooooooooooooooo 

000000000000000 000000 0000000000 00 000000000000000000 000000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 0000 0000000000000 
0000000 OOOO 00 OOOOOOOOOOOOOOOOOO 00000000000000000000 00000000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 000000 000000000000000 
000000000000000000000 00000000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 000000000000 000000000 

000000000 00 OOOOOOnOOOOOOOOOOOOOOOOOOOOO 0000 ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOO 00000000000000000000000000 OOOOOOOOOOOOOOOOOO 00 00 0000 0000000000000000000000000000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOvOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
000 00000000000000 00 0090 00 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOovOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00 OOOODOODOOOOO 
0000000 OOOOOOOOOOOCOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00 00000000 000000000 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 000000000000000000 00 00 ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00000000000000000000000000 OOOOOOOOOOOOOOOOOOOOOOOO 000000000000000000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 0000 0000 00 00 0000000000000 
000900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
000000009990000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
000 00 000000000000000000 000000 00 0000000000000000000000000000099000009CS9CC90000000999909999999999999999000900090000090000 
0000000000000090000000000000000000000 00000000000000000000 00 0090000000 00 00 000000 0000000000000000 000000 0000 00000000 0000000 
000000000 00 000000 00 00 00000000000000000000000900000000000000000000 OOOOOOOOOOOOOOOOOOOOOOOO 00000000000000 00000000000000000 
009909009000000900099090990000900000090000909000900000000000000000000000000000000000000000000000000000000000000000000000 
000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 090000000000000000000000000000000000 000000000000000000000000000 
0000000 0000000000000000000000 000000000000000000000000 OOOOOOOuwOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 000 
000 0000 0000000000000000000000 00000000000000000000 OOOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOO 00 0000000000000*^000000000000000000000 
000 OOOOOOOOOOOOOOOOOOOOOOOO 0000 00000000000000 000000090000000000000000000000000000000000000000000000000000000000000000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOoOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOO 000 
0000000 OOOOOOOOOOOOOOOOOO 0000000000000000 0000000000000000 00000000 OOOOOOOOOOOOOOOOOOOOOOOO 0000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000000000000090000000009000900000909000009900099000000000000000000000 
000000000000000090000000000000000000000000000000000000000000000000009000000000000000000000000000000000000000000000000000 
00090 00900900 90 000000000000000000000009000000900000000000000000000000000000000000000000000000000000000000000000000000000 
009000000000000000000000000000000000000000000000000000000000000000000QOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQ0090 
00000000090000000 990000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 0000 00000090 OOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 000000000000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOD 
000 OOOOOOOOOOt 0000000 00 OOOOOOOOOOOOOOOOOOOOOOOOOO 00 00 000000 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00 00 0000 90000000000 
00000000000 00000900000000000000 0000000090000000 0000009000 000000000000000000000000900000000000000000000000000000000000 000 
000 00 OOOOOOOCOOOOOOOOOOOO 000000000000 00000000000090000000000000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000090000090000009090990009000000909990000009009999099909999999999999999009990009 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooQoooooooooooo‘'oooooooooooooooooooooooooooooooooooooooooo90 
000000000 000900000000 0090 00 ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

990990099909000990099999909999090909909000000000000000000000000000000000000000000000000000000000000000000000000000000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00000 
00000000000000000 00 00 00 0000000000 oonooo 00 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 00 000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000009000000 00 9000 90000000000 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
OOoOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOC loooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
00000000000 00 00000000000000000000 000000000000000000000000000000000000000000000000000000000000000000000090000090000000000 
000000000000000000000000000000000000000000000000000000000000000000000000000000000000000009000999909990090900000990000000 
OOOOOOOOCOOOOOOOOOOOOOOOOOOOOOOOOOOCCCC 000000000000000000000000 0000000000000000000000000000000000000000000000 00000000000 
0000000 00 00 00 00000000000000000000000000000000000000000000000000000000000000000000000000000000000000900009000000000000000 
0000000000000000009999999990990900990990900999990000099000990900900000000000000000000000000000000000000000000DOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 0000 0000 0000000000 000 
000000000000000000000000000000000000000000000000000000000000000000009000000000000000000000000000000000000000000000000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO0OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOhOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
000000000000000000000000000000000000000000000000000000009000000000000000000000000000000000000000000000000000000000000000 
000000000000000000000000000000000000000000000000000000009099900099090099009999999990090000090000090Q09090090000900099900 
ooooooooooooooooooooooooooooooooooooooooouoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
00000 ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

09990990900000909099099D000099999009909000000D00000000000000000900000000000000000000000000000000000000000000000000000000 
00040000000000000C000000000000000000044000000000000000000000000900000000000000000000000000000000000000000000000000000000 
009400000000000000000000000000000000044000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0004D000000000000000000000000000000004400000000000000000000000000000000000000000000000000000000000000 0000000000000000000 
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This  is  a  copy  of  the  inputs  which  the  codes  prompts  the  user.  We  use  a  text 
file  which  was  piped  into  the  code. 

Stinp  ntmax-50 ,  dt  =  5  .  OOe-4 ,  cmax=0.0025,  eps  =  l.e-8, 
nout  =  l  Scend 
Scinp  models  &end 
£cinp  mode  =  0  £tend 
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D.  BOUK  INPUT 


Note  that  the  line  format  is  very  important  when  defining  the  problem 
geometry.  Comments  are  placed  within  the  input  file,  and  indicated  by  {*  ... 
*) .  However  the  comments  must  be  removed  prior  to  using  the  input  file. 

(*  General  definitions  and  gridding  for  the  problem  *) 

Big  RKA  Diode  /  Kyle  H  Setup  11/16/93  B  =  7.5  kG  E  =  400  kV 
RLIM  ZLIM  POTN  POT(POTN)  LSTPOT  MI  MAGSEG 

495  840  4  0.0,4.0E05,0.0,0.0  23-1 

(*  Axial  magnetic  field  for  each  axial  mesh  point  on  axis,  normalized  to 
MAGMLT.  The  values  begin  at  mesh  point  -6  and  goes  to  ZLIM+7 .  Notice  that 
1.0  is  used  to  denote  every  tenth  value.  This  helps  in  keeping  track  of  the 
magnetic  field  values.  The  data  are  actually  input  as  one  column,  however  it 
is  shown  here  in  multi-column  form.  *) 
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2.00 

0.99 

2 

224 

534 

0.99 

2.00 

2 

224 

600 

0.99 

2.00 

2 

224 

839 

0.99 

2.00 

2 

224 

840 

0.99 

0.00 

0 

223 

840 

2.00 

0.00 

0 

150 

840 

2.00 

0.00 

0 

0 

840 

0.00 

0.00 

0 

0 

800 

0.00 

2.00 

55 


0 


2 


0.00 


2.00 


0 

888 

UNIT  MAXRAY  STEP  NS  SPC  MASS  ZEND  VION  SAVE 
0.000332  27  0.1  13  0.6  0.0  840  -l.ES  0 

START  AV  AVR  RC  ZC  DENS  SURFACE  MAGMLT  MAGORD 
'GENERAL'  0  1.0  213.01  476.01  5000.0  1  7500.0  2 
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E.  RKA  INPUT 


The  Mathematica  calculations  for  the  space  charge  reduction  factor  are 
discussed  first.  The  equations  are: 


bsn/ :  bsn[(x_)  ,  (y_) ]  := 

(Pi*Sqrt [x*y] * (BesselY [0,  y] *BesselJ[0,  x]  -  BesselY[0,  x]*BesselJ[0,  y] ) ) / 
2 

Bsn/:  Bsn[(x_)  ,  {y_) ]  := 

(Pi*Sqrt [x*y] * {BesselY[l,  y] *BesselJ[l,  x]  -  BesselY[l,  x]*BesselJtl,  y]))/ 
2 

btn/:  btn[(x_)  ,  (y_) ]  := 

(BesseiY[0,  y]*BesselJ[0,  x]  -  BesselY(0,  x]*BesselJ[0,  y] ) / 

(BesselY[0,  y]*BesselJ[l,  x]  -  Besseiy[l,  x]*BesselJ[0,  y] ) 

btnh/ :  btnh[(x_)  ,  (y_) ]  := 

(BesselK[0,  x] *Besseil [ 0 ,  y]  -  BesselK[0,  y ] *BesselI [ 0 ,  x] ) / 

(BesselK[0,  y ] *BesselI [ 1 ,  x]  +  BesselKtl,  x] *BesselI [ 0 ,  y] ) 


g/: 


fl/: 

f2/: 

f3/: 

f4/: 


g[x_]  :=  x''(-2)  +  (Bessell[0,  c] /Bessell  [1,  c])*((btn[x  b  ,  x  c])''{-l)  - 
btnh[a  ,  b]/btn[x  c  ,  x  b] ) /x'^l  - 

(Bessell[0,  c]*Bsnix  b  ,  x  c]*btnh[a  ,  b] )/ (Bessell [1,  c]*bsn[x  b  ,  x  c] ) 
fl[x_]  :=  1  +  2*x*al  +  x''2*({l  +  al''2)  -  a2'-2)  +  2*x''3*al  +  x''4*al''2 

f2[x_]  :=  (1  +  2*x*al  +  x"2*((l  +  al^2)  -  a2''2))  -  2*x''3*al  +  x''4*al^2 

f3[x_]  :=  1  -  2*x*al  +  x''2*((l  +  al"'2)  -  a2^2)  +  2*x''3*al  +  x''4*al^2 

f4[x_]  :=  (1  -  2*x*al  +  x''2*((l  +  al^2)  -  a2''2))  -  2*x''3*al  +  x''4*al^2 


The  value  of  x  g[x]=  0  is  used  to  determine  the  space  charge  reduction 

factor.  The  problem  is  restricted  to  an  annular  electron  beam  in  a  hollow 
pipe.  The  values  of  the  radii  a  (wall  inner  radius),  b  (beam  outer  radius),  c 
(becim  inner  radius)  must  be  multiplied  by  Yo=(0/(Pc)  and  then  used  in  the 
equation  for  g. 

To  find  the  root  within  an  x-range  can  be  achieved  by: 

FindRoot [g [x_] , (x, xstart , xmin, xmax) ] . 

You  may  plot  g[x]  to  verify  that  there  is  a  root  by  using: 

Plot [g [x] , {x, xmin, xmax) ) . 

This  value  of  x  is  then  used  for  the  variable  a2,  while  al=fb/ ( fVy)  where  ft 
is  the  beam  plasma  frequency  and  f  is  the  rf  frequency.  The  roots  of  fl,  f2, 
f3,  and  f4  are  then  found  for  x  in  the  range  0  to  1 .  This  final  value  of  x  is 
the  plasma  reduction  factor  (PRF)  or  space-charge  reduction  factor  required  by 
RKA. 
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For  exeimple,  consider  the  case  of  a  300  kV,  10  kA  beeim  in  a  pipe  of  radius  a  = 
7.65  cm,  and  beam  position  b  =  6.6  cm  and  c  =  6.0  cm.  One  finds  the  following 
set  of  inputs  for  the  function  g[x] : 

7o  =  0.3441  cm'’  for  a  frequency  of  1.3  GHz 

a  =  2.6325 
b  =  2.2712 
c  =  2.0647 

g[x]=  0  yields  x=  2.2063  which  is  then  used  for  a2. 

Now  al  =  1.8775  yields  from  f4[x]=  0  that  the  PRF  =  0.248849.  This  is  due  to 
fl[x]  never  crossing  the  x-axis .  Also  f2  and  f3  are  not  physical  because  of 
the  crossed  +/-  signs.  That  is  the  signs  should  both  be  +  or  but  not  mixed. 


RKA  may  now  be  run  to  determine  the  beam  modulation  as  a  function  of  the  axial 
position  z.  A  seunple  deck  for  the  above  PRF  is  listed  below: 

Scparam 

ti=0.0,tf=1.0,nt=51, 

xi=0. 000005, xf =3 0. 000001, mx=151. 

d=l .56, f 0=1300. 0e+6,alpha=0. 050, 

volin=0 .300,i0=10.0e+3, rin=6 .0, rout =6 . 6, prf=0 .248849, 

nharm=2 , norbit=50 , 

nxint=5, 

yubhrm=l . 0 , nhrmint=4 , 
ylborb=0 . 0 ,yuborb=5 . 0, norbint=5, 
ylbvel=0 . 0 ,yubvel=l . 5, nvelint=5, 
y lbke=0 . 0 , yubke=2 . 5 , nkeint=5 , 
nequil=l, epslon=l .e-5, 
nrun=0 
/ 

One  finds  that  a  maximum  beam  modulation  of  10%  occurs  after  propagating  15 
cm,  13.44  cm  past  the  1.56-cm  modulating  gap.  Using  this  modulated  beam  as 
input  to  another  cavity  is  accomplished  by  renaming  the  generated  file 
RESTART. OUT  as  RESTART. IN,  changing  nrun=l,  changing  alpha  to  reflect  the 
voltage  induced  in  the  second  gap,  and  changing  xi  and  xf  for  the  positions  of 
the  second  gap  and  drift  space. 


58 


F.  KAGXC  INPUT 


Below  is  a  sample  input  deck  for  the  PIC  code  MAGIC.  The  plots  discussed  in 
section  3  can  be  found  in  the  input  listing. 

TITLE  •field(400kV  16kA)  cavity2cellshorter  ant8e6  “ ; 

COMMENT  "  Relativistic  Klystron  Amplifier  with  cup  cathode*; 

COMMENT  *  MAGIC  Jun  1993  SPARCIO  SIMULATION  FOR  KYLE  HENDRICKS  *; 

COMMENT  *  LAST  CHANGED  25  Jun  93  LESTER  BOWERS  (14Apr93)*; 
comment  "  /usrl2/home/les/magic/rka04/in04rka  incup*; 


DEFINE  PI  3.14159 

DEFINE  COURNT  . 7 
define  bi  16. E3 
define  bv  400. 
DEFINE  ANTRAMP  2.E-9 
define  dz  .001 


DEFINE  C  2.9979e+8;  DEFINE  bz  .6 

DEFINE  TSIM  59.e-9;  DEFINE  DUMP  29.E-9 

DEFINE  RAMP  lO.e-9;  define  delay  5.e-9 

DEFINE  AMPT  8.0e+6;  DEFINE  F  1.3E+9 

DEFINE  ROUND  l.e-7;  define  cavity  2 

define  dr  .001  ;  define  bkv=bv*1000. 


terminate  warning  ;  define  on  =  1  ;  define  off  =  0  ; 


define  idump  =  off 
define  iflux  =  off 
define  irange  =  off 
define  ivector  =  off 
define  icontour  =  off 
define  iantenna  =  on 
define  idisplay  =  on 
define  iobserve  =  on 
define  iphasespace  =  on 
define  iperspective  =  off 
define  ihighqfilter  =  off 
define  iemission_beam  =  off 


define  iemission_f ield  =  on 


define 

zO 

.0 

c 

start 

of 

rka 

define 

zl 

.1 

c 

end 

of 

shank 

define 

z2 

.2025 

c 

end 

of 

slant 

define 

z3 

.213 

c 

tip 

of 

cathode 

define 

z4 

.23 

c 

tip 

of 

anode 

define 

z5 

.267 

c 

end 

of 

a-k  gap 

define 

z6 

.377 

c 

start 

of 

first 

cavity 

define 

zl 

.387 

c 

start 

of 

first 

gap 

define 

z8 

.407 

c 

end 

of 

first 

cavity 

define 

z9 

.527 

c 

start 

of 

last 

cavity 

define 

zlO 

.537 

c 

start 

of 

last 

gap 

define 

zll 

.557 

c 

end 

of 

last 

cavity 

define 

zl2 

.66 

c 

end 

of 

rka 

define 

zm 

=  zl2 

define 

ngz 

12 

define 

rO 

.0 

c 

center  line 

define 

rl 

.07 

c 

cathode  tip 

define 

r2 

.074 

c 

cathode  tip 

define 

r3 

.076 

c 

anode  inside 

define 

r4 

.086 

c 

cavity  bottom 

define 

r5 

.  107 

c 

cathode  shank 

define 

r6 

.111 

c 

cathode  shank 

define 

r7 

.  107 

c 

cavity  top 

define 

r8 

.135 

c 

anode  break 

define 

r9 

.16 

c 

anode  outside 

define 

rm 

=  r9 

define 

ngr 

9 
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function  nxgrid{dxtotal , dxf irst , dxlast ) =2 . *dxtotal/ (dxf irst+dxlast ) ; 


define  zspectl  3.0  ;  define  zspect2  1.  ;  define  zspect3 

dzsuml  =  z4  ;  dzsum2  =  zll-z4  ;  dzsuin3  =  zm-zll  ; 


define  dzfirstl  =  zspectl*dz 
define  dzfirst2  =  dz 
define  dzfirst3  =  dz 


define  dzlastl 
define  dzlast2 
define  dzlast3 


define  ndzl  =  nxgrid (dzsuml , dzfirstl , dzlastl ) ; 
define  ndz2  =  nxgrid(dzsum2,dzfirst2,dzlast2) ; 
define  ndz3  =  nxgrid(dzsum3 , dzf irst3 , dzlast3 ) ; 
define  iO  2  ;  define  im  =  iO  +  ndzl  +  ndz2  +  ndz3 


dz 

dz 

zspect3*dz 


3.0 


define  aspectl  3.  ;  drsuml  =  rl  ; 

define  aspect2  1.  ;  drsum2  =  r7-rl  ; 

define  aspect3  3.  ;  drsum3  =  rm-r7  ; 

define  drfirst3  =  dr  ;  define  drlast3  =  aspect3*dr; 

define  drfirst2  =  dr  ;  define  drlast2  =  dr  ; 

define  drfirstl  =  aspectl*dr  ;  define  drlastl  =  dr  ; 

define  ndrl  =  nxgrid (drsuml, drfirstl, drlastl) ; 

define  ndr2  =  nxgrid(drsum2 , drf irst2, drlast2) ; 

define  ndr3  =  nxgrid(drsum3,drf irst3, drlastl ) ; 

define  kO  2  ;  define  km  =  kO  +  ndrl  +  ndr2  +  ndr3  ; 


DEFINE  BAREA  =  PI*(  R4**2  -  R2**2  ) 

define  bgamma  =  bv/511.  +  1. 

define  bbeta  =  sqrt (  1.  -  l./bgamma**2  ) 

DEFINE  BVEL  =  C*BBETA 

DEFINE  GAMV  =  BGAMMA*BVEL 

define  bj  =  bi/barea 

define  wavelen  =  100.*c/f 

define  rbeamavg=  .5*(rl+r2) 

define  twopir  =  2 . *pi*rbeamavg 

define  tempi  =  (bgamma** . 66667  -  l.)**l,5 

define  temp2  =  log (r3/rbeamavg) 

define  scli  =  8 . 55*templ/temp2 


c  kA  ; 


system  cylinder-theta; 
xlgrid  function  im  iO  zO 


ndzl  dzfirstl  dzsuml 
ndz2  dzfirst2  dzsum2 
ndz3  dzfirstl  dzsum3  ; 


x2grid  function  km  kO  rO  ndrl  drfirstl  drsuml 

ndr2  drfirst2  drsum2 
ndrl  drfirstl  drsuml  ; 


do  i  1, ngz, 1  ; 

engrid  xl  i'i'  z'i'  full  ; 
enddo  ; 

do  k  1 , ngr, 1  ; 

engrid  x2  k'k'  r'k'  full  ; 
enddo  ; 

z8_z6=z8-z6 ;  r7_r4=r7 -r4 ;  z8_z7=z8-z7 ;  r4_r3=r4-r3 ;  zll_z9=zll-z9 ; 

ungrid  xl  xli6  i6  full  ;  ungrid  xl  xli8  i8  full  ;  cavlnl=xli8-xli6  ; 
ungrid  xl  xlill  ill  full  ;  ungrid  xl  xli9  i9  full  ;  cavln2=xlill-xli9 
ungrid  x2  x2k4  k4  full  ;  ungrid  x2  x2k7  k7  full  ;  cavhit=x2k7-x2k4  ; 
ungrid  xl  xli7  i7  full  ,-  stubln=xli8-xli7  ; 
ungrid  x2  x2k3  k3  full  ;  stubht=x2k4-x2k3  ; 


conductor  cavities  anti-align 

i4  km  i5  k8  i5  k4  i4  k4  i4  k3 

i7  k3  i7  k4  i6  k4  i6  k7  i8  kl 
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i8  k3  ilO  k3  ilO  k4  i9  k4  i9  k7 

ill  k7  ill  k3  im  k3  im  km  iO  km  ; 

conductor  shank  align  iO  k6  il  k6  ; 

conductor  cathode  align  i3  kl  i3  k2  i2  k2  il  k6  ; 

conductor  noemit  align  i2  kO  i2  kl  i3  kl  ; 

define  dt  =  2 . *cournt/c*sqrt ( 1 . / (1 . /dz**2+l . /dr**2 ) )  ; 
define  fnyquist  =  l./(2.*dt)  ; 
define  ndt  =  tsim/dt  +  round  ; 
if  (  ihighqfilter  .eq.  on  )  then  ; 

define  cratio  =  4*cournt*cournt  ; 

z  define  gaa  =  .85  ; 

define  gaa  =  1.  -  1 . / (sqrt (3 . ) * (1 .+2 . *cournt) )  ; 

define  alphal  =  . 5*gaa/cratio  ; 

define  alpha2  =  1.  -  alphal  ; 

define  alpha3  =0.  ; 

define  iters  =  4  ; 

define  taul  =1.  ; 

define  tau2  =  .5  ; 

z  define  tau3  =  . 25* ( -3 . +sqrt (17 . ) )  ; 

z  define  tau4  =  . 25* { -3 . -sqrt { 17 . ) )  ; 

define  tau3  =  .2807764  ; 

define  tau4  =  -1.7807764  ; 

fields  all  biased  ndt  dt  alphal  alpha2  alpha3  iters  taul  tau2  tau3  tau4 
else  ; 

fields  all  biased  ndt  dt; 
endif  ; 

lookback  fields  all  1.  1.  anti-align  im  kO  im  k3  ; 

symmetry  axial  align  iO  kO  im  kO  ; 

define  dmaxnorm  l.e-4  ; 

define  nppc  2  ; 

define  mts  1  ; 

define  t2  =  delay  +  ramp  ; 

define  temp  .05  ; 

if  (  iemission_beam  .eq.  on  )  then  ; 

function  rho  data  4  0.0.  delay  0.  t2  bj  tsim  bj  ; 

emit  annular  bb  ; 

emission  annular  electron  nppc  mts 
beam  rho 

spacing  random  dmaxnorm  random 
thermal  gamv  temp  ; 

endif ; 

if  (  iemission_f ield  .eq.  on  )  then  ; 
define  threshold  1 . 5e7  ; 

'’"fine  breakdc'^m  1.  ; 

emit  annular  cathode  ; 

emission  annular  electron  nppc  mts 

field  breakdown  threshold  0.  0.  0. 
spacing  random  dmaxnorm  random 
energy  1 . e4  ; 
function  sph  =1.  ; 

function  voltramp  data  4  0.0.  delay  0.  t2  bkv  tsim  bkv  ; 

voltage  fields  tm  voltreimp  sph  1.  0.  1,  align  iO  k6  iO  km  ; 
endif ; 

CURRENTS  LCC  NO  NO  0 .  1.  ; 

FORCES  .5  1.  1.  ; 

KINEMATICS  ELECTRON  1  YES  NO  YES  EM  1  1  ; 

DEFINE  lk7  =  k7-l;  DEFINE  OMEGA  =  2.*PI*F  ; 
define  mzl  =  im-1;  define  k21  =  k2-l 

it  (  lantenna  .eq.  on  )  then  ; 
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FUNCTION  "TFUNIT)  =  MIN (T/ANTRAMP, 1 . ) *AMPT*SIN (OMEGA*T) "  ; 
antenna  tfun  null  null  el  16  18  lk7  lk7  ; 
endif ; 

c  bexternal  bz  .0  .0  ; 

preset  bird  read  MAGDATZ  surface  extmag  ascii  modify  replace  ; 
preset  b2rd  read  MAGDATR  surface  extmag  ascii  modify  replace  ; 

if  (  idump  .eq.  on  )  then 


DUMP 

FORMAT  ASCII  ; 

DUMP 

TYPE 

EOT  IND.ARY 

DUMP 

TYPE 

PHASES PACE 

DUMP 

TYPE 

PERSPECTIVE 

DUMP 

TYPE 

CONTOUR 

DUMP 

TYPE 

OBSERVE 

DUMP 

TYPE 

VECTOR 

DUMP 

TYPE 

RANGE 

DUMP 

TYPE 

FLUX 

DUMP 

TYPE 

GRID 

DUMP 
endif  ; 

TYPE 

PARTICLES 

DEFINE  FQN  .05E9  ;  DEFINE  FQX  5 . 9E9  ; 
output  meta  ; 
output  color  ; 


(  idisplay 

.eq. 

on  ) 

then; 

display  integer 

iO 

im 

kO 

km  nogrid  ; 

c 

full  device 

/ 

display 

real 

.36 

.43 

.04 

.13 

c 

first  cavity  ; 

display 

real 

.50 

.57 

.04 

.13 

c 

second  cavity  ; 

display 

real 

.19 

.29 

.05 

.15 

c 

a-k  gap  ; 

display 

real 

.09 

.22 

.02 

•15 

c 

cathode 

endif ; 

c  TAGGING  .8  ; 

DEFINE  NT  3  ;  DEFINE  FREQMAX  =  1 . / (2 . ♦NT*DT)  ; 

DEFINE  NSNAP  =  MAX (1 , DUMP / DT )  ;  DEFINE  NSTWO  =  2*NSNAP  ; 

TIMER  CHECK  PERIODIC  300  99999  600  ; 

TIMER  ENE-T  PERIODIC  0  99999  2  ; 

TIMER  SPIT  PERIODIC  NSNAP  99999  NSNAP  ; 

TIMER  FLX-T  DISCRETE  50  ; 

TIMER  RPUT  PERIODIC  NSTWO  99999  NSTWO  ; 

COURANT  SEARCH  ; 

DIAGNOSE  COURANT  101; 

DIAGNOSE  SPACING  100; 

DIAGNOSE  CONDUCTOR  101; 

STATISTICS  CHECK  ; 

if  (  irange  .eq.  on  )  then  ; 

c  across  bottom  of  cavity  gap  ; 

range  spit  1  field  el  i4  k5  i5  k5  ; 
range  spit  1  field  b2rd  iO  kl  im  kl  ; 

range  spit  1  field  bird  iO  kl  im  kl  ; 

endif  ; 

trajectory  99999  spit  1  electron  zO  zm  rO  rm  ; 

if  (icontour  .eq.  on)  then; 

contour  spit  field  el  iO  im  kO  km  boundary  yes  ; 

contour  spit  field  e2  iO  im  KO  km  boundary  yes  ; 

endif ; 

if  ( iperspective  .eq.  on)  then; 

perspective  rput  field  el  iO  im  kO  km  1  1  ; 
perspective  spit  field  bird  iO  im  kO  )<jn  1  1  ; 
perspective  spit  field  b2rd  iO  im  kO  km  11; 

endif ; 
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if  (iphasespace  .eq.  on)  then; 
phasespace  spit 

axes  xl  pi 

axis  X  .0  zm  .1 

axis  y  -.75e9  .75e9  ; 

endif ; 

if  (i vector  .eq.  on)  then; 
vector  spit  field  el  e2 
scale  log  1 
number  60  30 
axis  X  zO  zm  .1 
axis  y  rO  rm  .05  ; 

endif ; 

energy  ene-t  iO  im  kO  km  ; 

if  (iobserve  .eq.  on)  then; 
define  kf  5  ; 

c  observe  energy  voltage  1  0  interval  nt  fft  kf  ; 

c  window  frequency  fqn  fqx  ; 
observe  energy  electric  0  1  interval  nt  ; 

observe  energy  lookback  1  0  interval  nt  fft  kf 

window  frequency  fqn  fqx  ; 
c  E/M  field  energy  ; 
observe  energy  em  0  1  INTERVAL  NT  ; 
c  E/M  field  energy  rate  ; 

OBSERVE  ENERGY  em  0  0  INTERVAL  NT  fft  kf 
WINDOW  FREQUENCY  FQN  FQX  ; 

c  B_theta  one  cell  from  the  end  and  1  cell  below  Rbo  ; 
observe  field  b3  mzl  k21  mzl  k21  interval  nt  fft  kF 
window  frequency  fqn  fqx  ; 
c  voltages  -  first  gap  ; 

observe  field  el  i7  k3  i8  k3  interval  nt  fft  kf 
window  frequency  fqn  fqx  ; 
c  voltages  -  second  gap  ; 

observe  field  el  ilO  k3  ill  k3  interval  nt  fft  kf 
window  frequency  fqn  fqx  ; 

define  midgapl= ( i7+i8) /2  ;  define  midcavl= ( i8+i6 ) /2  ; 
c  Ez  field  in  the  middle  of  the  gap  on  top  of  the  beam  ; 
observe  field  el  midgapl  k2  midgapl  k2  interval  nt  fft  kf 
window  frequency  fqn  fqx  ; 
c  Ez  field  at  the  middle  botom  of  gap  ; 

observe  field  el  midgapl  k3  midgapl  k3  interval  nt  fft  kf 
window  frequency  fqn  fqx  ; 
c  Er  field  midcavity  one  cell  below  the  top  ; 
observe  field  e2  midcavl  lk7  midcavl  lk7  interval  nt  fft  kf 
window  frequency  fqn  fqx  ; 
c  voltage  on  exit  ; 

observe  field  e2  mzl  km  mzl  kO  interval  nt  fft  kf 
window  frequency  fqn  fqx  ; 

c  Ez  field  one  cell  below  the  outside  radius  of  the  beam  ; 
observe  field  el  mzl  k21  mzl  k21  INTERVAL  NT  FFT  KF 
window  frequency  fqn  fqx  ; 
c  injected  beam  current  ; 

observe  field  jl  i4  k3  i4  kO  interval  nt  ; 

c  from  first  gap  in  steps  of  3  cm  -  beam  current  ; 
define  stepl  .03  ;  define  istepl  =  stepl/dz  ; 
do  i  i8, im, istepl ; 
define  igl=  '  i  '  ; 

observe  field  jl  igl  k3  igl  kO  interval  nt  fft  kf 
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window  frequency  fqn  fqx  ; 

enddo ; 

if  (  iflux  .eg.  on  )  then; 

flux  zcavl  flx-t  all  indices  align  i6  kO  i6  k3  ; 
observe  flux  zcavl  current  ; 
flux  zslice  flx-t  all  plane  align  x  z7  ; 
observe  flux  zslice  current  ; 
endif ; 
endif ; 

START; 

STOP  ; 
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